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This  report  is  one  of  a group  of  Ship  Structure  Committee  Reports 
which  describes  the  SL-7  Instrumentation  Program.  This  program,  a 
jointly  funded  undertaking  of  Sea-Land  Service,  Inc.,  the  American 
Bureau  of  Shipping  and  the  Ship  Structure  Committee,  represents  an 
excellent  example  of  cooperation  between  private  industry,  regulatory 
authority  and  government.  The  goal  of  the  program  is  to  advance  under- 
standing of  the  performance  of  ships'  hull  structures  and  the  effective- 
ness of  the  analytical  and  experimental  methods  used  in  their  design. 

While  the  experiments  and  analyses  of  the  program  are  keyed  to  the  SL-7 
Containership  and  a considerable  body  of  the  data  developed  relates 
specifically  to  that  ship,  the  conclusions  of  the  program  will  be 
completely  general,  and  thus  applicable  to  any  surface  ship  structure. 

The  program  includes  measurement  of  hull  stresses,  accelerations 
and  environmental  and  operating  data  on  the  S.S.  Sea-Land  McLean, 
development  and  installation  of  a microwave  radar  wavemeter  for 
measuring  the  seaway  encountered  by  the  vessel,  a wave  tank  model  study 
and  a theoretical  hydrodynamic  analysis  which  relate  to  the  wave  induced 
loads,  a structural  model  study  and  a finite  element  structural  analysis 
which  relate  to  the  structural  response,  and  installation  of  long  term 
stress  recorders  on  each  of  the  eight  vessels  of  the  class.  In  addition, 
work  is  underway  to  develop  the  initial  correlations  of  the  results  of  the 
several  program  elements. 

Results  of  each  of  the  program  elements  are  being  made  available 
through  the  National  Technical  Information  Service,  each  identified  by 
an  SL-7  number  and  an  AD-  number.  A list  of  all  SL-7  reports  available 
to  date  is  included  in  the  back  of  this  report. 

This  report  contains  an  analysis  of  the  radar  wavemeter  and  the 
Tucker  wavemeter  systems  installed  on  the  Sea-Land  McLean.  Improve- 
ments in  the  systems  are  discussed  for  future  applications.  The  data 
presentation  procedure  to  be  used  in  subsequent  reports  is  discussed. 
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ABSTRACT 


So  that  more  precise  correlations  between  full  scale  observations 
and  analytical  and  model  results  could  be  carried  out,  one  of  the  objec- 
tives of  the  instrumentation  program  for  the  SL-7  class  container  ships 
was  the  provision  of  instrumental  measures  of  the  wave  environment. 

To  this  end,  two  wave  meter  systems  were  installed  on  the  S.S.  SEA-LAND 
McLEAN.  Raw  data  was  collected  from  both  systems  during  the  second 
(I973-197M  and  third  (197^-1975)  winter  data  collecting  seasons. 

it  was  the  purpose  of  the  present  work  to  reduce  this  raw  data, 
to  develop  and  implement  such  corrections  as  were  found  necessary  and 
feacible,  and  to  correlate  and  evaluate  the  final  results  from  the  two 
wave  meters.  In  carrying  out  this  work  it  was  necessary  to  at  least 
partly  reduce  several  other  channels  of  recorded  data,  so  that,  as  a 
by-product,  reduced  results  were  also  obtained  for  midship  bending 
stresses,  roll,  pitch,  and  two  components  of  acceleration  on  the  ship's 
br i dge. 

As  the  work  progressed  it  became  evident  that  the  volume  of  docu- 
mentation required  would  grow  beyond  the  usual  dimensions  of  a single 
technical  report.  For  this  reason  the  analyses,  the  methods,  the 
detailed  results,  discussions,  and  conclusions  are  contained  in  a series 
of  ten  related  reports. 

This  report  documents  some  background  analyses,  as  well  as  those 
which  were  necessary  to  develop  the  needed  corrections  to  the  raw  digi- 
tized data.  Implementation  of  the  results  of  the  various  analyses  in 
a second  data  reduction  stage  is  discussed. 


il 


CONTENTS 


page 

INTRODUCTION  I 

BACKGROUND  I 

PRELIMINARY  DISCUSSION  2 

Overall  Limitations  and  Goals  2 

ANALYSIS  OF  THE  RADAR  WAVE  HEIGHT  SYSTEM  4 

General  Analysis  4 

Axis  System 8 

Expansion  of  the  Slant  Range  Vector  11 

Expansion  of  Radar  Position  Vector 11 

The  Vertical  Component  of  Acceleration  14 

Relationship  Between  Ideal  Pendulum  Transducer 

Outputs  and  Euler  Angles  16 

Approximations  to  the  Required  Functions  of 

Roll  and  Pitch  Angles 18 

Error  Analysis:  Pendulum  Transducers  21 

Correction  of  Radar  Range  for  Vertical  Motion 

of  the  Ship . 29 

The  Low  Frequency  Cut-Off 31 

Filtering  Methods  - Sampled  Data  36 

Preliminary  Numerical  Development; 

Fast  Convolution  Method  38 

Radar  Wave  Processing  System 39 

Simulation  of  Data  for  Final  Development  of 

Low  Frequency  Cutoff  Parameters  of  H(u'e) 41 

Final  Version  of  the  Low  Frequency  Cutoff  Parameters  42 

Results  of  Final  Precision  Checks  44 

Sample  Evaluation  with  Actual  Data  47 

Retrospective  52 

ANALYSIS  OF  THE  TUCKER  METER 54 

Background 5^ 

The  SL-7  Tucker  Meter  Installation  56 

Analysis  of  Possible  Corrections  57 

A Study  of  Errors  Induced  by  Double  Integration  57 


i i i 


page 

INDIRECT  METHODS  OF  WAVE  MEASUREMENT 67 

THE  FINAL  DATA  REDUCTION  AND  PRESENTATION  PROCEDURE  69 

REFERENCES 75 


LIST  OF  FIGURES 


PAGE 

1 OVERALL  GEOMETRY  - RADAR  6 

2 SEMI-FIXED  AND  MOVING  AXIS  SYSTEMS  (NO  HEAVE  10 

OR  SWAY) 

3 RADAR  AIMING  NOTATION  12 

4 RELATIVE  POSITIONS  OF  ACCELEROMETERS  AND  RADAR  12 

5a  SAMPLE  TIME  HISTORY,  RUN  401  26 

5b  SAMPLE  TIME  HISTORY,  RUN  401  26 

6 LOCATION  OF  INTEGRATION  ERRORS  IN  THE  WAVE  30 

FREQUENCY/HEADING  PLANE,  SHIP  SPEED  20  KT. 

7 LOCATION  QF  INTEGRATION  ERRORS  IN  THE  WAVE  30 

FREQUENCY/HEADING  PLANE,  SHIP  SPEED  30  KT. 

8 LOCATION  OF  INTEGRATION  ERRORS  FOR  QUARTERING/  30 

FOLLOWING  SEAS  IN  THE  WAVE  LENGTH  - HEADING  PLANE 

9 DEFINITION  OF  H(we)  39 

10  PRECISION  OF  DOUBLE  INTEGRATION  BY  FAST  CONVOLUTION  41 

METHODS  (HARMONIC  AND  NARROW  BAND  ACCELERATION) 

11  PRECISION  OF  DOUBLE  INTEGRATION  BY  FAST  CONVOLUTION  41 

METHODS  (WIDE  BAND  ACCELERATION) 

12a  SAMPLE  RANGE,  ACCELERATION,  DISPLACEMENT,  AND  43 

ENCOUNTERED  WAVE  TIME  HISTORIES,  RUN  401 

12b  SAMPLE  RANGE,  ACCELERATION,  DISPLACEMENT,  AND  43 

ENCOUNTERED  WAVE  TIME  HISTORIES,  RUN  401 

13  RUN  401,  RADAR  WAVE  SPECTRUM  AND  COMPONENTS  45 

14  COMPARISON  OF  NOMINAL  MIDSHIP  MOMENT  RESPONSE  FROM  47 

SPECTRAL  RATIO  WITH  MODEL  DATA  (RUN  401) 

15  DIAGRAM  OF  THE  TUCKER  METER  INSTALLATION  53 

16a  RUN  401,  SAMPLE  TIME  HISTORIES,  TUCKER  METER  AND  58 

MEAN  DYNAMIC  HEAD 

16b  RUN  401,  SAMPLE  TIME  HISTORIES,  TUCKER  METER  AND  58 

MEAN  DYNAMIC  HEAD 

17  SPECTRA  DERIVED  FROM  TUCKER  DATA  (RUN  401)  58 

18  RUN  401  — VOYAGE  32W  — TAPE  145  — INDEX  17  — INTERVAL  1 63 

19  RUN  401  --  VOYAGE  32W  — TAOE  145  — INDEX  17  — INTERVAL  1 64 


SHIP  STRUCTURE  COMMITTEE 


The  SHiP  STRUCTURE  COMMITTEE  is  constituted  to  prosecute  a research 
program  to  improve  the  hull  structures  of  ships  by  an  extension  of  knowledge 
pertaining  to  design,  materials  and  methods  of  fabrication. 


RAD.M  W.  H.  Banker t,  USCG  (Chairman) 
Chief,  Office  of  Merchant  Marine  Safety 
U.S.  Coast  Guard  Headquarters 


Mr.  P.  M.  Palermo 
Asst,  for  Structures 
Naval  Ship  Engineering  Center 
Naval  Ship  Systems  Command 

Mr.  John  L.  Foley 

Vice  President 

American  Bureau  of  Shipping 


Mr.  M.  Pitkin 
Asst.  Administrator  for 
Commercial  Development 
Maritime  Administration 

Mr.  C.  J.  Whitestone 
Engineer  Officer 
Military  Sealift  Command 


SHIP  STRUCTURE  SUBCOMMITTEE 

The  SHIP  STRUCTURE  SUBCOMMITTEE  acts  for  the  Ship  Structure  Committee 
on  technical  matters  by  providing  technical  coordination  for  the  determination 
of  goals  and  objectives  of  the  program,  and  by  evaluating  and  interpreting  the 
results  in  terms  of  ship  structural  design,  construction  and  operation. 


NAVAL  SEA  SYSTEMS  COMMAND 

Mr.  R.  Johnson  - Member 

Mr.  J.  B.  O'Brien  - Contract  Administrator 

Mr.  C.  Pobler  - Member 

Mr.  G.  Sorkin  - Member 

U.S.  COAST  GUARD 

LCDR  T.  H.  Robinson  - Secretary 
LCDR  S.  H.  Davis  - Member 
CAPT  C.  B.  Glass  - Member 
Dr.  W.  C.  Dietz  - Member 

MARITIME  ADMINISTRATION 

Mr.  F.  Dashnaw  - Member 
Mr.  N.  Hammer  - Member 
Mr.  R.  K.  Kiss  - Member 
Mr.  F.  Seibold  - Member 

MILITARY  SEALIFT  COMMAND 

Hr.  T.  W.  Chapman  - Member 
COR  J.  L.  Simmons  - Member 
Mr.  A.  B.  Stavovy  - Member 
Mr.  D.  Stein  - Member 

AMERICAN  BUREAU  OF  SHIPPING 


NATIONAL  ACADEMY  OF  SCIENCES 
SHIP  RESEARCH  COMMITTEE 

Mr.  0.  H.  Oakley  - Liaison 
Mr.  R.  W.  Ruroke  - Liaison 

SOCIETY  OF  NAVAL  ARCHITECTS  fi 
MARINE  ENGINEERS 

Mr.  A.  B.  Stavovy  - Liaison 
WELDING  RESEARCH  COUNCIL 
Mr.  K.  H.  Koopman  - Liaison 

INTERNATIONAL  SHIP  STRUCTURES 
CONGRESS 

Prof.  J.  H-  Evans  - Liaison 
U.S.  COAST  GUARD  ACADEMY 
CAPT  W.  C.  Nolan  - Liaison 

STATE  UNIV.  OF  N.Y.  MARITIME  COLLEGE 
Dr.  W.  R.  Porter  - Liaison 
AMERICAN  IRON  & STEEL  INSTITUTE 
Mr.  R.  H.  Sterne  - Liaison 


Mr.  S.  G.  Stlansen  - Chairman 
Or.  H.  Y.  Jan  - Member 
Mr.  I.  L.  Stern  - Member 


U.S.  NAVAL  ACADEMY 
Dr.  R.  Bhattacharyya  - Liaison 

U.S.  MERCHANT  MARINE  ACADEMY 
Or.  Chin-Bea  Kim  - Liaison 


vi 


INTRODUCTION 


In  the  analysis  of  the  wave-induced  ship  hull  strain  data  obtained 
by  SSC  in  the  1960's  it  was  necessary  to  infer  the  wave  environment  from 
estimated  Beaufort  wind  speeds.  An  extraordinary  amount  of  work  was 
required  to  develop  the  inferential  techniques.  These  techniques  appear 
to  suffice  for  valid  prediction  of  long-term  trends  because  a great  deal 
of  averaging  is  carried  out.  Unfortunately  when  verification  of  short- 
term statistical  predictions  is  desired,  the  use  of  wind  as  a wave 
environment  index  appears  to  be  less  than  satisfactory. 

As  a consequence  it  was  one  of  the  objectives  of  the  SL-7  full- 
scale  instrumentation  program  to  provide  a direct  instrumental  measure 
of  the  wave  environment  so  that  more  precise  correlations  could  be  made 
between  full-scale  observations,  and  analytical  and  model  results. 

To  this  end  the  ship  was  fitted  with  a micro-wave  radar  relative 
wave  meter  and  various  motion  sensing  devices.  A "Tucker  Meter"  pressure 
actuated  wave  height  sensing  system  was  also  installed. 

The  purpose  of  the  present  project  is  to  reduce  and  analyze  the 
resulting  wave  meter  data  obtained  on  the  SEA-LAND  McLEAN  in  the  second 
(1973-1974)  and  third  (1974-1975)  winter  recording  seasons. 

The  purpose  of  the  present  report  is  to  document  the  various 
analyses  carried  out  so  as  to  produce  the  basic  data  reduction  system, 
and  to  describe  that  system. 


BACKGROUND 


In  the  documentation  of  the  present  project  it  has  been  necessary 
for  practical  reasons  to  assume  on  the  part  of  the  reader  a general 
familiarity  with  the  Ship  Structure  Committee's  SL-7  measurement  program. 
The  primary  background  references  for  the  present  project  are  References  1 
through  4.  Reference  1 is  the  basic  documentation  of  the  full-scale 
instrumentation  system.  References  2 and  3 contain,  for  both  recording 
seasons  in  question,  a quite  full  account  of  instrumentation,  basic  record- 
ing, and  the  nominal  circumstances  surrounding  the  present  data.  These 
references  also  contain  results  of  analyses  of  longitudinal  vertical  mid- 
ship bending  stress  which  were  carried  out  according  to  the  methods  of 
Reference  4. 

The  first  report  under  the  present  project  is  Reference  5.  This 
reference  contains  the  detail  of  the  selection  of  the  data  set  for  the 
1973'1974  winter  season,  the  methods  utilized  in  the  digitization,  the 
calibration  of  the  data,  and  results  of  a first  analysis  of  the  data. 


PRELIMINARY  DISCUSSION 


Overall  Limitations  and  Goals 

The  present  state-of-art  of  prediction  or  extrapolation  of  ship 
motions  or  stresses  is  grounded  firmly  in  the  linear  short-term  random 
approach  as  first  put  together  by  St.  Denis  and  Pierson  (References  6,7). 
Under  this  model  the  seaway  is  characterized  by  a directional  variance 
spectrum,  and  is  assumed  to  be  statistically  stationary  and  Gaussian. 
Because  the  ship  responses  are  assumed  to  be  linear  functions  of  the 
seaway,  they  are  also  stationary  and  Gaussian,  and  may  be  characterized 
by  variance  spectra. 

As  a practical  matter,  the  ultimate  results  of  any  full-scale 
instrumentation  program  are  interpreted  in  light  of  the  linear  random 
model.  In  the  data  reduction  of  ship  responses  the  absolute  minimum 
requirement  is  to  produce  the  statistical  variance,  or  its  square  root 
(rms),  or  a theoretically  related  statistic  such  as  the  average  or  root- 
mean-square  deviation  of  the  maxima  of  a record  from  zero.  When  simple 
time  histories  are  involved  as  raw  material,  the  data  reduction  cake  is 
often  usefully  iced  with  various  statistics  of  the  maxima  and  minima  of 
the  record  (the  extremes,  the  average  of  10th  highest,  etc.).  Finally, 
for  some  correlation  purposes,  estimates  of  the  variance  spectra  may  be 
desired  and  can  be  obtained  with  standard  techniques. 

In  general,  when  the  measurement  of  sea  waves  from  a moving  ship 
is  considered,  some  considerable  complication  enters  the  data  reduction 
picture.  The  complication  arises  from  two  sources: 

a.  Short  crestedness  of  the  seaway. 

b.  Speed  of  advance  of  the  ship. 

In  the  case  of  the  wave  measuring  instrumentation  installed  in  the 
SEA-LAND  McLEAN,  the  best  that  can  be  expected  is  that  the  reduced  data 
will  represent  t;ie  wave  elevations  at  a point  moving  with  the  ship.  Thus, 
even  at  best,  the  two  complications  just  mentioned  are  of  concern  in  the 
present  problem. 


According  to  the  basic  model,  the  directional  spectrum  of  encounter 
is  equal  to  the  directional  wave  spectrum  times  a frequency  mapping.  The 
frequency  mapping  simply  converts  the  point  of  reference  from  a fixed 
point  (wave  frequency)  to  a moving  point  (encounter  frequency).  Because 
the  waves  do  not  change  with  arbitrary  changes  in  coordinate  systems, 
the  frequency  mapping  preserves  variance.  According  to  the  assumptions 
in  the  basic  model,  the  integration  of  the  directional  spectrum  over 
direction  and  frequency  is  the  variance.  Furthermore,  an  adequate  esti- 
mate of  the  variance  may  be  obtained  with  knowledge  of  the  time  variation 
of  wave  elevation  at  a single  point.  Accordingly,  there  is  no  basic  theo- 
retical obstacle  to  the  estimation  of  the  seaway  variance  from  elevations 
at  a movi ng  point. 
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Within  the  f-ame  of  reference  of  the  basic  statistical  model  the 
same  arguments  hola  for  averages  of  maxima  of  wave  elevation  time  histo- 
ries, with  the  interpretative  exception  that  the  resulting  estimates  of 
wave  elevations  (or  heights)  must  be  considered  to  be  indicative  only  of 
the  rise  and  fall  along  the  track  of  the  ship  — no  information  about 
short  or  long  crestedness  is  involved  and  no  information  about  wave  direc- 
t i ona 1 i ty . 

Thus  the  time  history  of  wave  elevation  at  a moving  point  can  theo- 
retically yield  reasonable  estimates  of  wave  variance  or  average  maxima. 
There  are  some  circumstances  where  this  is  not  practically  so.  Consider 
the  situation  in  following  or  quartering  seas  where  the  predominant  wave 
energy  is  concentrated  in  a narrow  frequency  and  directional  band  and 
the  typical  component  of  wave  celerity  in  the  direction  of  ship  advance 
is  equal  to  ship  speed  (a  500-ft,  10-sec  following  wave  for  30-kt  ship 
speed,  for  instance).  In  this  circumstance,  the  apparent  wave  elevation 
at  the  moving  point  will  be  slowly  varying  or  practically  constant  for  a 
short  sample,  and  the  typical  variance  estimate  (which  involves  correc- 
tion to  zero  sampl e mean)  can  be  appreciably  in  error.  This  is  a sampling 
problem  rather  than  theoretical,  but  is  nonetheless  real  and  must  be  con- 
sidered. There  will  undoubtedly  be  some  "loss"  in  variance  in  following 
sea  samples  regardless  of  the  distribution  of  variance  in  the  actual 
directional  wave  spectrum. 

There  is  a much  higher  order  of  complication  in  the  spectrum  of 
wave  elevations  at  a moving  point.  Supposing  a true  time  history  to  be 
available  from  a wavemeter  --  the  spectrum  of  this  function  of  time  takes 
the  following  form: 

U(o)  ) = S.  (u>  ,X  )dX 
v e J 1 v e’  e e 

+ I Ss(vXe)dXe 

+ J S3(VXe)dXe  0) 

where 

U(u)e)  = spectrum  of  encounter  at  a moving  point 

S.(u)g,Xe)  = directional  spectrum  of  encounter,  mapping  region  i 

U)e  = encounter  frequency 

Xg  = wave-ship  heading 

To  map  the  actual  directional  wave  spectrum  into  the  encounter 
plane,  it  must  be  divided  into  three  regions  of  uu  and  X where  unique  one- 
to-one  mappings  exist.  In  general,  there  may  be  a contribution  to  the 
spectrum  at  a given  encounter  frequency  from  three  different  wave  fre- 
quencies. As  shown,  the  directional  spectrum  is  imbedded  in  an  integral 
equation.  Even  if  only  one  region  is  involved  in  the  actual  directional 
spectrum  (a  fairly  narrow  directional  spread  in  head  seas  for  instance) 
there  will  be  insufficient  data  contained  in  the  time  history  of  one  of 
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the  wavemeters  under  discussion  to  derive  any  directional  properties.  The 
scalar  wave  spectrum  (integration  over  heading)  is  derivable  if  the  ship 
speed  is  zero.  In  summary,  the  directional  wave  spectrum  is  a function  of 
time  and  space  (frequency  and  direction).  Multiple  wavemeters  are  required 
at  zero  ship  speed  to  provide  even  the  possibility  of  obtaining  directional 
spread.  The  scalar  wave  spectrum  is  not,  in  general,  derivable  with  only 
one  piece  of  data  if  the  ship  is  at  speed. 

Clearly,  it  will  not  be  possible  to  extract  everything  that  could 
be  desired  from  the  present  wavemeter  data  — a limitation  which  was 
accepted  at  the  very  outset  of  the  program.  However,  the  overall  full- 
scale/theoretical/model  correlation  problem  is  such  that  instrumental 
measurement  of  even  the  simplest  wave  statistic  represents  a potentially 
significant  improvement,  and  thus  the  extraction  from  the  data  of  these 
simple  statistics  is  the  basic  minimum  goal  of  the  present  project. 


ANALYSIS  OF  THE  RADAR  WAVE  HEIGHT  SYSTEM 


General  Analysis 

The  characteristics  of  the  radar  system  (the  "OWHS  Radar")  are 
detailed  in  Reference  8.  Essentially,  the  radar  yields  the  range  from 
the  radar  dish  to  the  water  surface.  Because  the  radar  is  aimed  at  an 
angle  with  respect  to  ship  co-ordinates,  the  radar  output  is  called  the 
"slant  range."  The  radar  system  carried  out  only  one  further  step  in 
processing,  and  that  is  to  subtract  a nominal  "initial"  range  from  the 
instantaneous  range  before  carrying  out  a digital  to  analog  conversion. 

The  initial  nominal  range  upon  ship's  departure  was  not  recorded.  This 
was  apparently  because  the  system  has  automatic  features  which  insure 
automatic  acquisition  of  the  initial  range  — and  re-acquisition  in  case 
of  temporary  return  signal  loss. 

Thus  the  radar  yields  only  a measure  of  the  relat ive  motion  between 
ship  and  wave.  The  main  effort  in  the  present  radar  data  reduction  process 
is  to  correct  this  relative  range  for  angular  and  vertical  motion  of  the 
ship,  and  in  so  doing  arrive  at  an  estimate  of  the  encountered  wave  ele- 
vat i ons . 

Conceptually,  the  correction  of  the  radar  range  for  ship  angular 
orientation  is  a straight-forward  vector  operation.  Unfortunately,  the 
nature  of  the  radar  system  makes  it  impossible  to  establish  the  length 
of  the  range  vector  solely  from  the  instrumental  data.  The  only  reason- 
able option  appeared  to  be  to  assume  that:  a)  the  sample  mean  of  the 
radar  data  in  each  interval  corresponds  to  the  position  of  the  nominal 
still  water,  zero  speed,  waterline  of  the  ship;  b)  to  correct  the  instru- 
mental data  to  this  mean;  and  c)  to  add  to  the  corrected  data  the  distance 
from  radar  antenna  to  still  water  as  computed  from  the  ship's  departure 
drafts.  If  the  sample  mean  is  really  reasonably  close  to  the  still  water 
waterline  this  procedure  might  result  in  total  range  errors  of  2 -*  4 feet 
out  of  75.  Since  the  corrections  are  of  a cosine  nature  it  is  thought 
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that  the  final  results  will  not  be  too  far  off  the  mark.  Unfortunately, 
however,  it  is  not  possible  to  find  out  just  how  far  wrong  the  above 
procedure  might  be  in  any  given  case. 

As  noted  in  Reference  5,  the  signal  re-acquisition  feature  also 
apparently  produced  large  sudden  shifts  in  the  mean  value  of  the  radar 
signal.  It  was  thought  best  to  try  to  eliminate  instances  of  this 
behavior  from  the  data  set  rather  than  attempt  analysis. 

To  simplify  the  analysis  it  was  assumed  that  the  range  had  been 
corrected  for  the  initial  bias  as  noted  above.  The  situation  may  then 
be  outlined  in  a general  way  as  in  Figure  1,  in  which  the  static  (in 
harbor)  case  is  to  the  left  and  the  dynamic  case  is  shown  to  the  right. 

There  are  three  types  of  hardware  in  the  system,  all  fixed  in  the 

ship. 

1 . The  radar 

2.  An  accelerometer  package 

3.  Angular  motions  transducers. 

The  accelerometer  package  is  intended  to  yield  the  vertical  translation 
of  the  ship  at  that  point.  It  is  thus  helpful  to  establish  a fixed 
vertical  reference  plane  parallel  to  the  mean  water  surface  (normal  to 
local  gravity)  which  passes  through  the  accelerometer  package  when  the 
ship  is  at  rest  in  calm  water.  This  situation  is  indicated  at  the  left 
of  Figure  1.  The  reference  plane  is  (H-T)  above  mean  water  level,  where 

H = Height  of  accelerometer  package  above  baseline 
T = Local  ship  draft. 

In  the  first  four  voyages  in  the  second  season  the  radar  and  the 
nearest  accelerometers  are  not  located  at  the  same  point  on  the  ship. 

This  fact  is  expressed  in  the  left  side  of  Figure  1 by  the  vector  (J 
which  is  defined  to  be  a fixed,  three-dimensional  position  vector  for 
the  radar  relative  to  the  accelerometer  package. 

Note  that  it  has  been  convenient  to  define  positive  vertical 
directions  down . For  this  reason  and  for  convenience  in  analysis  the 
position  of  the  radar  package  is  shown  above  the  radar  though  this  is 
not  generally  so  in  the  ship. 

The  radar  is  aimed  in  a three-dimensional  sense  and  so  the  slant 
range  for  the  static  case  is  a constant  position  vector  S. 

Adopting  the  notation: 

of  Q 

downward) 
of  S 

downward) 


= vertical  component 
(pos i t i ve  d i rected 

Sz  = vertical  component 
(pos i t ive  d i rected 
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CAL  REFERe'OGE’  PLAver 


FIGURE  I OVERALL  GEOMETRY  - RADAR 


Thus : 


(2) 


(h-t)  = gz  + s2 

(by  def i n i t i on) 

Now  proceeding  to  the  dynamic  case  (right  side  of  Fig.l),  the  ship 
is  assumed  to  translate  and  rotate.  The  reference  planes  are  as  before, 
but  the  local  water  surface  deviates  from  mean.  The  wave  elevation,  a 
function  of  time,  is  defined  as  £(t) , positive  upward  (to  correspond  to 
the  usual  ideas  about  the  ups  and  downs  of  waves) . 

The  quantity  2^(0  is  the  vert i ca I translation  of  the  accele- 
rometer package. 

The  position  of  the  radar  relative  to  the  accelerometer  is  expressed 
by  the  vector  P(t) . This  vector  is  the  same  as  ij  relative  to  the  ship,  but 
is  in  a different  direction  relative  to  fixed  co-ordinates  due  to  angular 
motion  of  the  ship. 

Finally  the  slant  range  is  defined  to  be  ft(t) . Relative  to  ship 
co-ordinates  the  slant  range  has  the  same  d i rect i on  as  ?.  Relative  to 
fixed  coordinates  the  direction  is  different  according  to  the  ship's  angu- 
lar mot i on . 

The  vertical  components  of  the  vectors  are: 

Pz(0  = vertical  component  of  position  vector  -- 
pos i t ive  down. 

R (t)  = vertical  component  of  slant  range  vector  — 
positive  down. 

The  component  Pz(t)  is  a function  of  time  through  the  time  varia- 
tion of  pitch  and  roll.  The  component  Rz ( t ) is  a function  of  time  for 
the  same  reasons  --  in  addition  to  the  variation  of  1 R(t)  ! , which  is  the 
quantity  measured. 

All  of  the  parts  of  the  problem  have  now  been  indicated  in  a gen- 
eral way.  From  Figure  1 the  general  expression  for  the  wave  elevation 
may  be  written: 

-C(t)  = Zt(t)  + Pz(t)  + Rz(t)  - - Sz  (3) 

Equation  (3)  is  the  basic  framework  for  the  data  reduction  process. 
Because  the  geometry  is  three-dimensional  the  expressions  are  deceptively 
simple  with  respect  to  implementation.  The  additional  analysis  required 
may  be  divided  into  several  problem  areas  of  which  the  first  few  appear 
as  fol lows : 

I.  Provision  of  consistent  axis  systems  and  the  required 
transformations  between  ship  and  fixed  co-ordinates. 
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2.  Expansion  of  vectors  P(t)  and  R(t)  into  forms  compatible 
with  the  data,  and  computation  of  ( t)  and  Rz(t) 


3.  Analyses  of  angle  transducers. 

4.  Consideration  of  the  generation  of  the  vertical  component 
of  acceleration  and  its  double  integration. 

Axis  System 

To  start,  an  axis  system  fixed  in  the  ship  and  one  fixed  in  space 

may  be  adopted  in  accordance  with  the  Eulerian  conventions  in  common  use, 

Reference  9.  It  is  supposed  that  a right  handed  set  of  axes  is  fixed 

relative  to  the  earth  o that  the  X , Y axes  are  horizontal  and  the  Z 

° ° o 

axis  is  directed  downward.  The  set  of  axes  fixed  in  the  ship  (x,y,z)  are 
so  chosen  that  it  is  initially  coincident  with  the  fixed  axes. 

The  x axis  is  positive  forward,  the  y axis  is  positive  to  starboard, 
and  the  z axis  is  positive  downward. 

In  order  to  express  the  position  of  a point  in  the  ship  in  terms  of 
the  fixed  co-ordinate  system  define: 

= (X  ) co-ordinate  of  position  of  origin  of  ship  axes 

Yt  = (Y  ) co-ordinate  of  position  of  origin  of  ship  axes 

= (Zq)  co-ordinate  of  position  of  origin  of  ship  axes  (4) 

and  (x1,y1,z1)  = co-ordinates  of  the  point  with  respect 
to  the  ship  axes 

Then  the  co-ordinates  of  the  point  in  terms  of  the  fixed  axis  system  are: 


X!  = Xt 

+ 

an*i 

+ 

+ 

ai3  Z1 

Y:  = Yt 

+ 

aslxi 

+ 

+ 

aS3Z! 

Zt  =Zt 

+ 

am*i 

+ 

a3Syi 

+ 

a3  3Zl 

where  the  a.,  are  the  direction  cosines  between  the  ship  axes 
and  the  fixecl  axes. 

Following  the  cited  conventions,  the  direction  cosines  are  derived 
from  three  operations: 

1.  Rotate  the  X , Y axes  about  Z an  angle  of  yaw  *. 

o'  o o 3 

2.  Rotate  the  result  about  the  body  axis  y through 
an  ang 1 e of  pitch,  9 . 

3.  Rotate  this  result  about  the  body  axis  x through 
the  angle  of  rol  1,  cp. 
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The  right  hand  rule  is  followed  so  that  yaw  is  positive  as  the  bow 
moves  starboard,  pitch  is  positive  as  bow  moves  up,  and  roll  is  positive 
starboard  side  down.  The  direction  cosines  resulting  from  these  conven- 
tions are: 

DIRECTION  COSINES  (a..)  OF  BODY  AXES  RELATIVE  TO  FIXED 
AXES  IN  TERMS  of  0,  9,  t|r. 


J-a 

J-3 

i-1 

COS0  COSl|l 

-coscp  sim|i  + sin0  sincp  cosi|i 

sincp  sin\|f  + sin©  coscp  cosijr 

i=2 

cos©  sin'll 

cos  9 cositr  + sin0  sin9  sim|t 

-sincp  cos'll  + sin©  coscp  sin^i 

i=3 

-s  i n0 

cos  9 s i n9 

cos©  cos 9 

Now  in  the  present  problem  the  interest  is  in  the  vertical  compo- 
nents of  motion,  and  in  any  event  there  is  not  enough  data  being  acquired 
to  make  yaw  a useful  input.  Accordingly,  the  earth-fixed  co-ordinate 

system  X , Y , Z is  of  little  direct  use  in  the  problem  and  a semi-fixed 
o o o 

coordinate  system  to  provide  vertical  reference  is  more  to  the  point. 

If  the  ship  is  assumed  always  to  travel  in  the  direction  of  positive  X 
and  the  translation  X{  of  Eq . (4)  is  assumed  to  be  ship  velocity  times 
time,  the  relations  between  the  new  semi-fixed  axis  system  and  the  ship 
axes  are  found  by  substituting  X^  = 0 and  = 0 in  Eq.  (5)  and  in  the 
above  table  of  direction  cosines. 

Defining  co-ordinates  of  a point  in  the  semi-fixed  system  as 
X,  Y,  Z the  result  for  a point  (x,y,z)  in  the  ship  system  becomes: 


X 

= 

bnx 

+ 

bi*y 

+ 

•v. 

Y 

= Yt 

+ bsix 

+ 

bssy 

♦ 

bS3 

Z 

“ Zt 

+ b3ix 

+ 

b^y 

+ 

b33 

and  the  direction  cosines  b..  are  as  follows: 

U 

DIRECTION  COSINES  (b..)  OF  BODY  AXES 
U 

RELATIVE  TO  SEMI-FIXED  AXES  IN  TERMS  OF  0,  9 


j-1 

j*2 

j-3 

i = l 

cos© 

sin©  s i n9 

sin©  cos9 

i=2 

0 

cos  9 

-s  in9 

i-3 

-s  i n0 

cos  © s i n9 

cos©  cos 9 

Figure  2 indicates  the  sign  conventions  and  the  order  of  intro- 
ducing angular  motion  (pitch-roll)  for  zero  translation  of  the  ship  axis 
system. 
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Expansion  of  the  Slant  Range  Vector 


The  orientation  of  the  radar  antenna  is  specified  by  an  azimuth 
relative  to  the  bow  and  a depression  angle  from  the  horizontal. 


Referring  to  Figure  3,  the  azimuth  angle  (y)  is  defined  as  posi- 
tive for  a swing  to  starboard.  The  depression  angle  (6)  will  be  con- 
sidered positive  in  the  direction  of  the  z axis.  Defining  direction 
cosines  of  the  vector  ? (and  1?  relative  to  ship  coordinates)  with  respect 
to  ship  co-ordinates  x,y,z  as  cx , c^ , cg : 


Then 

ci 

= 

cosy  cos 6 

c 

= 

siny  cos 6 

2 

C3 

= 

s i n6 

which  is  to  say 

: 

S 

X 

= 

1 * 1 - | 

?! 

cosy 

cos  6 

S 

= 

me  = i 

?! 

s i nV 

cos6 

y 

2 

s 

= 

] ? ! c3  = ' 

? 1 

s i n6 

z 

3 

Now  defining: 

i “*  i 

S 

= 

1 s i 

Ra(t)  = | R(t)  | 

(S  and  R (t)  correspond  to  the  absolute  range  from  the  radar) 

O 3 


(7) 


(8) 

(9) 


Following  Eq . (6)  and  considering  only  the  relations  between  the  antenna 
and  the  end  of  the  ranges: 


S c, 
o 3 


Rz  (t)  = Ra(t)|-c1  sin8  + cg  cos9sincp  + cg  cosScoscpj. 
= 8a(0[cab3i  + c b + c3b33] 


(10) 


(ID 


Thus  the  absolute  range  of  the  radar  is  separated  from  the  functions  of 
ship  roll  and  pitch. 


Expansion  of  Radar  Position  Vector 

It  will  be  most  convenient  to  specify  the  location  of  the  radar 
relative  to  the  accelerometer  in  terms  of  heights  and  offsets.  Accord- 
ingly the  components  of  the  position  vector  (J  are  outlined  in  Figure  4. 
Spec i f i cal  1 y, 

Q = distance  of  radar  forward  of  accelerometers. 

(positive  forward) 
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“y 


athwartship  distance  of  radar  from  accelerometer, 
(positive  to  starboard) 


= vertical  distance  of  radar  below  accelerometers 
positive  for  down  (below),  negative  for  radar 
above  accelerometers. 

Then  from  Eq.  (6)  as  before: 


Pz(t)  = b^  + b3s(iy  + 

= -Q^sinG  + Q^cosQsincp  + Q^cosGcos1? 


The  Vertical  Component  of  Acceleration 


(12) 


The  accelerometer  packages  installed  in  the  ship  are  bi-axial 
(vert i ca 1 -athwartsh i p)  and  each  accelerometer  is  of  the  seismic  type. 


Seismic  accelerometers  are  sensitive  to  angular  orientation  of 
their  sensitive  axis  in  that  a signal  component  due  to  gravitational 
attraction  is  produced  which  is  generally  indistinguishable  from  true 
acceleration.  In  particular,  the  apparent  acceleration  (in  g's)  due  to 
gravity  is  equal  to  the  cosine  of  the  angle  between  the  accelerometer's 
sensitive  axis  and  true  vertical. 

For  analysis  purposes  a tri-axial  accelerometer  package  will  be 
assumed.  The  sensitive  axes  of  the  accelerometers  are  assumed  aligned 
parallel  to  the  ship  fixed  (x,y,z)  axes  (Fig. 2).  The  signals  from  the 
three  accelerometers  are  defined: 


A (t) 
Ay  (t) 
Az(t) 


= output  from  accelerometer  aligned  with  x-axis 
= output  from  accelerometer  aligned  with  y-axis 
= output  from  accelerometer  aligned  with  z-axis 


Positive  acceleration  will  be  defined  as  being  in  the  positive 
direction  of  the  x,y,z  axes  and  the  dimensions  of  all  outputs  are  assumed 
to  be  in  g's. 


It  is  conventional  to  electronically  offset  the  output  of  a 
vertically  oriented  accelerometer  to  compensate  for  constant  gravity  — 
that  is  to  consider  the  output  as  zero  when  the  ship  has  no  heel  and  is 
at  rest  in  calm  water.  Accordingly,  for  present  purposes  the  apparent 
acceleration  signal  from  the  accelerometer  on  the  z axis  will  be 

(Az(t)  - 1) 

The  cosines  of  the  angles  between  the  sensitive  axes  of  the  accel- 
erometers and  the  ZQ  axis  (vertical)  are  given  in  the  table  following 
Eq.  (6).  Thus  the  component  of  the  sensed  accelerations  due  to  gravity 
are  (in  g's): 
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along  x 

axis: 

-b3  1 

= +s i n9 

y 

axis: 

-b3S 

= -cos9sin9 

z 

axis : 

-b33 

= -cos9coscp 

(Note  that  the  negative  sign  is  because  deflection  of  accelerometer  mass 
is  opposite  to  that  defined  as  positive  acceleration.)  Thus  the  sensed 
accelerations  corrected  for  the  influence  of  gravity  are: 


x axis  — A^(t)  + b3i 
y axis  — Ay(t)  + b3jj 
2 axis  — Az(t)  - l+b33 

When  there  is  no  dynamic  component  and  no  list  or  trim,  the  total  z-com- 
ponent  = -I  = gravity  component;  = -b3g  ; and  A^ ( t)  = 0.  The  above  three 
components  define  a corrected  translational  acceleration  vector  defined 
in  ship  co-ordinates.  For  purposes  of  double  integration,  the  component 
along  true  vertical  is  required  and  this  is  obtained  by  exercising  Eq . (6): 

= b3TAx(t)  + b3SAy(t)  + b33AzM  ~ (b93-’)  f 3) 

In  the  present  instrumentation  setup  there  are  no  x-axis  accel- 
erometers. What  is  assumed  is  that  the  mass  of  the  ship  is  so  large 
relative  to  longitudinal  excitation  forces  that  the  true  longitudinal 
acceleration  is  zero.  Exercising  Eq.  (6)  and  simplifying,  the  true 
longitudinal  acceleration  (component  on  x axis)  becomes: 

AL(t)  = bn\(t)  + blzAy(t)  + b13Az(t)  ' bta  (,4) 

Assuming  A^(t)  = 0 and  manipulating  Eq . (14)  an  estimate  for  the  un- 
measured component  becomes: 

A (t)  = - -g— — A (t)  ~(A  (t)  - 1)  (15) 

* l l Y n 

Substituting  in  Eq . (13)  there  results  the  following  estimate  of  true 
vertical  acceleration: 


Av(t)  = Ay(t)/b3s  - b3l|  + Az(t){b33 


1 -b. 


33 


and  upon  substitution  of  the  expressions 


for  b.  . 
•J 


and  some  reduction: 


A (t) 
v 


— — ^ A (t)  + 
cos 0 y' 


coscp 
cos  9 


Az(t) 


^~coscp 

l_cos9 


'] 


(16) 


Equation  (16)  is  then  the  final  result  for  the  vertical  accelera- 
tion under  the  assumption  that  no  true  longitudinal  acceleration  is 
present . 
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Regardless  of  the  assumption  about  true  longitudinal  acceleration 
the  true  transverse  acceleration  (A^t))  turns  out  to  be: 

A^Ct)  = costPAy(t)  - s incp(Az(t)  - 1)  (17) 

Relationship  Between  Ideal  Pendulum  Transducer  Outputs  and  Euler  Angles 

In  the  present  case  the  pitch  and  roll  transducers  are  of  the 
pendulum  type  and  are  fixed  in  the  ship.  In  each,  the  pendulum  is  con- 
strained to  move  in  a plane  normal  to  its  axis.  For  pitch  the  plane  of 
motion  is  the  x-z  plane,  (rotation  about  y axis)  and  for  roll,  rotation 
about  the  x axis  in  the  y-z  plane.  If  the  pendulum  transducer  is  working 
ideally  (or  statically)  the  pendulum  aligns  itself  with  the  intersection 
of  its  plane  of  freedom  and  the  plane  in  which  lie  both  its  rotation  axis 

and  the  semi-fixed  vertical,  Z . 

o 

Referring  to  Figure  2,  and  assuming  the  roll  pendulum,  the  pendulum 
is  free  to  rotate  in  the  y-z  plane  about  the  x axis  and  is  in  equilibrium 
along  the  intersection  of  this  plane  with  the  x-Z  plane.  This  intersec- 
tion is  noted  z1  in  Figure  2 by  virtue  of  the  convention  about  pitching 
the  ship  first.  The  direction  cosines  of  this  line  with  respect  to  the 

fixed  X , Y , Z axes  are 
o'  o'  o 

(sin0,  0,  cosG) 

The  angle  indicated  by  the  pendulum  is  that  between  z and  z1  in 
Figure  2 since  in  the  zero  or  up-right  position  of  the  ship  9=0  and 
z and  z1  are  coincident.  The  direction  cosines  of  z are: 

(sin9  cos 9,  -sin9,  cos9cos9) 

Accordingly  the  cosine  of  the  angle  sensed  by  the  pendulum  is: 

sin29cos9  + cos2 9cos9  = cos9 

and  the  ideal  roll  pendulum  output  is  the  same  as  9 in  the  previous  sec- 
tion, and  the  sign  convention  will  be  assumed  to  be  the  same. 

Now  considering  pitch,  the  plane  of  rotation  of  the  pendulum  is 
defined  by  y=0.  The  direction  cosines  of  the  normal  to  the  plane 
(the  y axis)  are: 


(sin9sin9,  cos9,  cos9sin9) 

Accordingly,  the  direction  cosines  of  the  pendulum  (£,m,n)  at  any  position 
in  its  plane  of  rotation  must  satisfy: 

0 = ■£  sin9sin9+  m cos9  + n cos9sin9 

At  equilibrium  the  pitch  pendulum  must  also  lie  in  the  y-Z  plane. 
This  plane  is  not  coincident  with  the  Yq-Zo  plane  because  of  the  roll 
angle  and  thus  the  ideal  pitch  pendulum  output  is  not  the  same  as  the 
Euler  pitch  angle  9. 


15 


Assume  the  direction  cosines  of  the  normal  to  the  y-Z  plane  are 

o 

(p,q,r).  The  normal  to  the  plane  must  also  be  normal  to  the  Zq  axis 
which  results  in  r=0.  The  normal  to  the  plane  is  also  normal  to  the 
y axis  and  this  requirement  means: 

0 = p sin9sin9+  q cos9  + r cos 9s  in'? 

Finally  the  direction  cosines  must  satisfy  the  relationship 

s s , s , 
p + q + r =1 

With  these  relations  the  direction  cosines  of  the  normal  to  the  y-Z  plane 
become: 


P = 'Z., 


+ sin?9tan2<T 


q = -sin9tancp/y  1 + sin29tanS<p 


From  this  result  it  follows  that  the  direction  cosines  of  the  equilibrium 
position  of  the  pendulum  (Z,m,n)  must  satisfy 

0 = £-m  sin9tan<p 
They  must  also  satisfy 

.s  s a , 

1 + m + n =1 

as  well  as  the  equation  derived  previously  for  the  x-z  plane: 

0 = 1 sin9sin<P  + m cos^P  + n cos9sin<P 
Solution  of  the  three  relations  simultaneously  yields 

„ -s i n9cos9tan9s inT 
1+s  i ns  9tanScp 

-cos 9s  i ncp 


+J  1 +s  i n3  9tanScp 
n = coscp ^ 1+s  i nS  9tan?cp 

The  angle  sensed  by  the  pitch  pendulum  will  be  called  (o)  and  it 
is  the  angle  between  the  equilibrium  position  defined  above  and  the  z axis. 
Accord i ngly, 

coso  = £b13  + mbg3  + nb33 


= cos  e/y  1+s  i nS  0tanScp 


Manipulating  Eq.  (18),  there  results  the  relations  between  the  observed 
angle  a and  the  requi  red  angle  9 (for  1 co  ! < n/2)  : 
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cos  9 


s i n0 


cosa 


cost^/ ( 1 + cos3  Ot  an2  cp) 


s i no 


+/  l+cos3otan3  cp 


09) 


Clearly  for  zero  or  very  small  roll  o and  9 are  nearly  the  same, 
which  is  according  to  reason. 


The  preceding  development  allows  the  functions  of  the  Euler  angles 
to  be  expressed  in  terms  of  the  ideal  observable  angles. 

It  is  assumed  that  the  sense  of  the  pendulum  output  is  the  same  as 
noted  in  Figure  2: 

cp  = Roll  angle  as  observed  by  pendulum  transducer  — 
positive  for  starboard  side  down. 

o = Pitch  angle  as  observed  by  pendulum  transducer  -- 
positive  for  ship  bow  up. 

The  functions  required  by  the  development  in  previous  sections  are: 


cr 

GJ 

II 

. „ -s  i no 

-sin0  = — ■.  ■ ■ 

+ / 1+cos  atan3^ 

(20) 

^3  2 

. . „ cosotancp 

cosAsintp  = — — — - 

+/  l+cos2otan‘ cp 

(21) 

^3  3 = 

cosor 

COS  GcOSCp  = ; 

/ 2 2 
{/  1+cos  otan  cp 

(22) 

/S  i nca 
lcose' 

sincpcos9  /,.  s..  s 

= + / 1+cos  atan  cp 

COSO  v/ 

(23) 

(COSCP) 

vcos9' 

cos?<P  / , , s 

= / 1+cos  otan  cp 

COSO  J 

(24) 

Approximations  to  the  Required  Functions  of  Roll  and  Pitch  Angles 

Equations  (20)  through  ( 24 ) are  likely  to  be  unnecessarily  time 
consuming  to  evaluate  and  a reasonable  simplification  is  desirable.  The 
pitch  and  roll  transducers  have  stops  at  ± 45°  so  that  neither  angle 
sensed  can  be  greater  under  any  circumstances. 

According  to  the  typical  results  indicated  in  Reference  5,  a "large" 
osciMatory  pitch  is  ± 2 or  3 degrees  and  a "large"  oscillatory  roll  is 
± 15  . Considering  natural  limits,  the  magnitude  of  the  largest  pitch 
observed  in  model  tests  at  D.L.  (under  unusually  severe  conditions)  is 
± 15  . In  the  present  case  a 10°  pitch  moves  the  bow  up  (or  down)  about 
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90  feet  or  3 times  draft.  Clearly,  an  approximation  good  to  ± 10°  pitch 
seems  more  than  adequate  in  the  present  case.  Similar  natural  limits 
upon  roll  are  not  so  obvious.  Plus-or-minus  15°  is  out  of  the  usual  small 
angle  approximation  range  and  was  observed  in  the  fourth  voyage.  Most 
seamen  are  willing  to  recount  tales  of  90  out-to-out  rolls  (±45  degrees) 
It  appears  conservative  to  split  the  difference,  and  so  the  objective  as 
far  as  roll  was  concerned  was  to  make  an  approximation  good  to  ± 30°,  and 
at  least  acceptable  to  ± 45°. 

Toward  this  end  the  expressions  in  o and  cp,  Eqs.  (20)  through  (24), 
were  expanded  in  series  retaining  terms  to  5th  or  6th  order.  Finally,  in 
each  case  the  three  or  four  terms  of  lowest  order  were  retained.  These 
approximations  were  compared  with  the  exact  evaluations  of  Eqs.  (20)  to 
(24)  over  a ± 10  range  of  o and  a ± 45°  range  of  9.  Some  minor  adjust- 
ment of  coefficients  was  made  to  improve  the  correspondence  with  the  exact 
evaluations  for  large  angles  and  the  final  approximations  are  as  follows: 


-s  i no 


b3  1 

\ 

j l+cosSQ'tanScp 

h - - 

cosotancp 

.3  2 

\ 

/ 2 2 
/ l+cos  otan  cp 

b - - 

COSO 

3 3 

> 

/ 2 2 

J l+cos  otan  cp 

s i nc f 
cos  9 

s i n9coscp  / , 
/ l+cos 

COSO  v 

-o  + O.I6670  + 0.5ocpS  (25) 

cp  - 0.16679  - 0. 439oS  (26) 

1 _ 0.499^  - 0.5os  + 0.7  oV  (27) 


cp  - 0.166793  + 0.43cpo® 


coscp  cos  cp  r s s 

+ = — / 1+cos  otan  cp 

cos  9 cosO  yj 


1 - 0.49^  + 0.5oa  - 0.7oscp2 


(28) 


(29) 


The  range  of  the  possible  values  of  the  coefficients  in  Eqs.  (25) 
to  (29)  is  in  the  interval  ±1  . The  maximum  absolute  error  in  the  approx- 
imations indicated  is  ± 0.001  over  the  ranges  -10  ^ o £ 10°,  -30  s ¥ s 30° 
The  absolute  error  enlarges  to  ± .003  to  ± 0.006  up  to  roll  angles  of  40° 
and  is  ± .004  to  ± 0 . 0 1 3 at  45°  roll. 


The  specified  resolution  of  the  radar  is  roughly  1 foot.  If  the 
approx imat ions  above  are  reasonable  they  should  introduce  errors  much 
smaller  than  1 foot  for  perfectly  measured  data.  The  magnitude  of  the 
range,  Rg(t),  is  in  the  neighborhood  of  75  feet.  If  depression  angles 
(6)  between  60  and  70  are  chosen  with  azimuths  (y)  between  45°  and  90°, 
the  coefficients  ci , c^ , c3 , Eq.  (7),  vary  roughly  as  follows: 
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.87  < c3  < .94 


In  the  range:  roll  = ±30°  and  pitch  = ±10°  the  magnitudes  of  b , b32,  b33 
are  .16,  .7  and  1.  at  most  (with  the  worst  combinations  of  signs),  to  that 
if  each  of  the  coefficients  bgl,  b32,  b33  is  in  error  by  +0.001,  the 
maximum  error  in  the  multiplier  of  R (t)  comes  out  to  be  about  0.0015  and 
this  typically  corresponds  to  an  error  of  0.1  feet.  Errors  in  the  approx- 
imations are  not  usually  of  the  same  sign  so  that  the  expected  maximum 
error  in  the  corrected  range  is  significantly  less  than  this.  When  roll 
angles  are  greater  than  30  the  corresponding  error  is  significantly 
greater,  but  the  waves  themselves  would  be  expected  to  be  large  and  the 
incluence  of  the  approximation  error  relatively  smaller. 

Pretty  much  the  same  argument  applies  to  the  position  component, 
(Pz(t)-Q.  ) Eq.  (12).  In  this  case  the  possible  use  of  the  midship  accel- 
erometer package  is  excluded  because  it  has  not  been  recorded  on  the  same 
tape  as  the  radar,  roll  and  pitch.  In  any  event  it  would  be  a poor  choice 
since  Q,  would  be  =«  500  feet  and  an  error  of  1/2  in  9 would  result  in 
a 5 foot  error  in  the  position  component.  Of  the  two  accelerometer 
packages  remaining  (Refs.  1,2,3)  the  one  farthest  from  the  radar  has: 

Q*  * 47  ft 

% - 50  ft 

d «=  -40  ft 
z 

and : 

Pz(t)-Qz  * 47  b31  + 50  b3S  - 40  (b33-l) 

The  coefficient  (b  -1)  is  always  negative  and  at  most  is  about  -.15. 

The  coefficients  b3,  and  b3g  are  at  most  .16  and  .7  so  that  for  errors 

of  +.001  in  each  coefficient  the  maximum  error  in  (P  (t)-Q  ) due  to  the 

z z 

approximations  appears  to  be  about  0.05  feet.  In  the  actual  data  reduc- 
tion case  the  accelerometer  package  nearest  the  radar  was  used  and  thus 
the  potential  error  is  even  less  than  that  just  cited. 

Turning  to  the  acceleration  A^(t),  Eq.  ( 1 6) , when  roll  is  small 
the  vertical  accelerometer  dominates,  its  coefficient  is  unity  and  the 
absolute  .001  approximation  error  translates  to  0.1%.  When  roll  is  large 
the  horizontal  acceleration  has  been  observed  to  dominate  and  in  this  case 
the  error  might  range  between  0.2  and  0.5%. 

It  should  be  emphasized  that  all  the  above  argument  involve*  per- 
fectly measured  data.  It  appears  that  errors  or  possible  lack  of 
resolution  in  the  basic  range,  acceleration  and  angle  data  may  have  a 
much  greater  influence  upon  error  in  the  final  product  than  the  errors 
inherent  in  the  approximations,  Eqs.  (25)  to  (29). 
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Error  Analysis:  Pendulum  Transducers 


The  pendulum  transducers  are  specified  to  have  a 2 Hz  or  higher 
natural  frequency  with  what  appears  to  be  between  0.4  and  0.7  times 
critical  damping  depending  upon  ambient  temperature.  Accordingly  the 
errors  due  to  the  response  of  the  instruments  themselves  should  be  neg- 
ligible below  0.2  Hz  and  at  least  acceptably  small  below  0.4  Hz.  It 
appeared  from  References  2 and  3 that  little  if  any  roll  and  pitch 
response  occurs  at  frequencies  greater  than  0.3  Hz;  that  is,  that  the 
ship  itself  does  not  respond  as  a rigid  body  to  excitations  by  waves  which 
are  short  enough  to  result  in  encounter  frequencies  above  0.3  Hz.  Because 
of  the  simplicity  of  the  instrument,  correction  of  data  for  the  pendulum 
frequency  characteristics  would  be  feasible,  but  it  appeared  reasonable 
to  ignore  the  pendulum  dynamics  altogether. 

The  above  considerations  do  not  remove  all  the  potential  dynamic 
problems  from  the  pendulum  transducers.  For  frequencies  sufficiently 
below  pendulum  resonance  the  pendulum  aligns  itself  along  a line  opposite 
to  the  direction  of  the  resultant  translational  force  on  its  pivot.  If 
the  pivot  is  not  accelerating,  this  line  is  the  local  vertical.  All  the 
previous  derivations  of  ideal  pendulum  angles  have  been  for  the  static 
case  or  for  the  case  of  no  translational  acceleration  of  the  pivot.  The 
purpose  of  the  present  development  is  to  assess  the  errors  introduced  by 
translations  of  the  ship  in  way  of  the  pendulum  transducers. 

First  it  will  be  assumed  as  in  a previous  development  that  the 
inertia  of  the  ship  is  so  large  in  relation  to  longitudinal  excitation 
that  fore  and  aft  accelerations  of  the  ship  at  the  pendulum  location  are 
zero.  It  is  known  that  the  ship  heaves  and  highly  likely  that  it  sways 
so  that  two  translational  accelerations  may  be  assumed: 

A^.(t)  = Horizontal  acceleration  in  the  direction  of 
the  (fixed)  Yq  axis.  Figure  2.  (in  g's) 

A (t)  = Vertical  acceleration  in  the  direction  of 
the  (fixed)  Z^  axis.  Figure  2.  (in  g's) 

(in  the  following  the  functional  notation  will  be  dropped  and  the  time 
variation  will  be  understood.) 

With  these  assumptions  and  sense  conventions,  the  Y component  of 
the  t rans 1 at i ona 1 force  on  the  pendulum  pivots  is: 

(w  Ar) 

the  component  is  zero, 
and  the  Zq  component  is: 

w (A  -1) 
v 

(where  w is  the  pendulum  weight) 
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Thus  the  apparent  gravity  which  the  pendulum  "sees"  is  a vector  with: 


Xq  component 


Yq  component  = -AT 
Zq  component  = 1 -A 

and  which  has  direction  cosines 


p = 0 

0 ■ -s/J * "A/ 

r - ^ * ('-V2 

As  in  the  previous  development  the  location  of  the  pendulum  is  along 
the  intersection  of  a plane  perpendicular  to  its  axis  and  a plane  passing 
through  the  axis  and  the  local  gravity  vector  above.  Once  the  direction 
cosines  of  the  pendulum  position  (.C,m,n)  are  found  the  cosine  of  the  angle 
the  pendulum  makes  with  the  z axis  can  be  found  as  in  a previous  section. 
This  process  was  carried  out  for  the  pitch  and  roll  pendulums  and  the 
results  are  as  follows: 

For  the  roll  pendulum: 


COS0 


tancp 


cos(3  = cosCf 


(30) 


For  the  pitch  pendulum: 


cost  = 


cos9  - tan<P 

r 


/’ + + s i r?  9)  tanScp  - 2 cosGtancp 


where: 


(30 


P = the  angle  sensed  by  the  roll  pendulum 
t = the  angle  sensed  by  the  pitch  pendulum 

Substitution  of  Eq . ( H)  and  some  purely  algebraic  manipulation  serves 
to  put  Eqs.  (30)  and  (31)  in  the  following  form: 


tanp 


tancp  + A 
I - Atancp 


(32) 


tanT 


tana- 

1 - Atancp 


(33) 
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where 


A = a £211 


p /T7 


? 2 
cos  atari  9 


coso 


r cos  0 


Also  following  from  Eqs . (32),  (33): 


tan(P-cp) 

tan(T-o) 


A tancpcosas  i no 
1 - A tancpcos'o 


(34) 

(35) 

(36) 


These  last  two  expressions  are  for  the  errors  (t-cv,P-9)  which  are 
introduced  by  the  translational  accelerations  of  the  pendulum.  As  written, 
the  R.H.  sides  of  Eqs.  (32)  to  (36)  involve  onl\  ideal  parameters  — none 
of  which  can  be  observed  directly  by  the  installed  instrumentation. 


Ideally,  there  is  the  possibility  of  eliminating  the  unkown  ideal 
parameters  by  using  the  outputs  from  the  accelerometer  package  mounted 
alongside  the  pendulums.  Equations  for  the  true  vertical  and  transverse 
accelerations  are  Eqs.  ( 1 6)  and  (17).  Combining  these: 


2 

r 

and  thus 

A 


With: 


-At  A - (Az-I)tan9 

1 -A  = COSP  (A  -1)  + A tan9 
v z y 


H - tan<y 
1 + Htan9 
A 


Substitution  of  Eq . (38)  into  Eq . (32)  yields: 


(37) 


(38) 


A 

tanp  = H = (39) 

z 

No  expression  for  the  unknown  angle,  9,  results  because  both  the  accel- 
erometers and  the  pendulum  are  seismic  instruments.  Both  "see"  the  same 
effective  local  gravity  vector  and  are  both  recording  the  same  information. 
To  correct  sensed  pendulum  angles  to  true  it  is  necessary  to  have  an 
independent  measure  of  the  true  transverse  acceleration,  To  correct 
sensed  accelerations  to  true  vertical  and  transverse  acceleration  it  is 
necessary  to  have  an  independent  measure  of  true  angles. 

With  the  instrumentation  installed  there  is  no  way  to  make  the 
corrections  to  indicated  angles  which  are  implied  by  Eqs.  (32)  to  (36). 
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As  a practical  matter  the  pitch  errors  shown  in  Eq.  (36)  are 
probably  fairly  small.  Making  the  small  angle  approximation: 


T-cr 


Acp 


A true  transverse  acceleration  of  the  ship's  C.G.  of  O.lg  might  be  quite 
a lot,  and  this  coupled  with  10°  roll  yields  errors  of  the  order  of  2%. 

If  roll  and  sway  are  roughly  in  quadrature  as  expected,  the  error  will  be 
smaller.  It  should  be  remarked  that  the  potential  first  order  errors  in 
pitch  angle  were  assumed  away  at  the  outset  by  assuming  A^(t)  = 0 (Eq.l4). 

The  error  situation  with  roll  is  not  pleasant.  A true  transverse 
acceleration  of  O.lg  should  correspond  roughly  to  5 error,  which  could 
easily  be  relatively  very  large. 

With  the  data  at  hand  it  must  be  assumed  that  actual  sway  accel- 
erations are  negligible  in  way  of  the  pendulums  --  there  is  no  way  to 
check  this  assumption  with  the  data  which  has  been  acquired. 

Recapitulation  of  Analysis  Thus  Far  and  Sample  Evaluation 

Substituting  Eqs.  (II)  and  (12)  in  Eq . (3),  the  expanded  expression 
for  the  wave  elevation  becomes: 

-C(t)  - \(0  + (P  (t)-Q  ) + (A.(t)-S  ) 

v z z z z 

= Z.(t)  + b„d_  ♦ b,aa..  + (b„-l)a. 


't 
+ R 


3*  X 


“S3 


,<*>  h 


b„,  + c b„  _ + 

3 1 s 3 S 


(M» 


The  translation  Zt(t)  is  the  double  integration  of  the  true  vertical 
acceleration  at  the  accelerometer  package  location,  Eq.  (16): 


V'v  [V’] 


Hnf  A (0  + A (t)  - \ssq.  r 

cosa  v ' cosy  1 ' 1 cos  a 


cose  "y"'  cos9  2 Lcos9  J 

The  notation  is  summarized  for  convenience  as  follows: 


(M) 


A.  Parameters  Fixed  by  Installation  or  Ship  Geometry  (Numerical 
values  were  noted  in  Ref.  5 for  the  second  season). 

Q.  — The  coordinates  of  the  radar  position  relative 

’ to  the  accelerometer  position. 

S — The  vertical  position  of  radar  antenna  above 

waterl i ne. 

c. ,c„,c.  --  Direction  cosines  of  radar  beam  relative  to 

1 7 2 7 3 

ship-fixed  co-ordinate  system. 
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B.  Recorded  Functions  of  Time 


Ra(0 

Az(t),Ay(t) 

b.  . 

•J 


The  indicated  slant  range,  the  actual  range, 
not  the  variation  of  the  range  about  some  mean 
value. 

Vertical  and  horizontal  acceleration  signals 
as  recorded  by  the  fixed  bi -axial  accelerom- 
eter package. 

Functions  of  roll  and  pitch  angles  cp  and  9; 
approximations,  Eqs . (25)  through  (27). 


The  first  four  terms  of  Eq,  (40)  together  are  the  variation  in 
vertical  position  of  the  radar  antenna  from  its  nominal  height  above 
water,  S^.  The  b..  coefficients  are  complicated  functions  of  time  depen- 
dent angles.  The  slant  range,  P (t),  is  the  entire  range  — not  the 
deviations  from  some  mean. 


Because  Eqs.  (40)  and  (41)  involve  products  of  measured,  time 
dependent,  variables  a great  part  of  the  implementation  had  to  be  in  the 
time  domain.  In  particular,  there  was  no  reasonable  alternative  to  the 
computation  of  the  quantities: 


C-C(t)  - zt(t)]  = Rc(t) 

and;  A (t) 

v 

in  the  time  domain. 

The  "wave  elevation"  in  Eq.  (40)  is  not  necessarily  zero  mean. 

The  development  thus  far  has  been  in  terms  of  coordinates  fixed  for  zero 
ship  speed.  When  the  ship  is  underway  there  will  be  small  trims  and  heels 
constant  for  sample  durations  and  probably  some  bodily  sinkage.  In  addi- 
tion there  may  also  be  some  dynamic  swell-up  of  water  at  the  radar's  point 
of  aim.  These  things  add  up  to  a constant  bias  in  the  vertical  position 
of  the  radar  and  a bias  in  the  slant  range  itself.  There  is  in  addition 
the  problem  cited  at  the  beginning  of  the  analysis  of  the  radar  system; 
that  is,  the  absolute  value  of  the  slant  range  cannot  be  known  with  great 
accuracy.  The  assumptions  for  coping  with  this  were  outlined  previously, 
and  make  no  provision  for  sinkage,  dynamic  swell-up,  etc.  Long  term 
apparent  biases  cannot  reasonably  be  handled  in  the  integration  for  Z(t) 
because  of  the  discontinuous  nature  of  the  records.  However,  a zero  bias 
in  the  derived  wave  elevation  is  not  serious  since  all  the  useful  results 
from  the  final  product  involve  an  explicit  or  implicit  correction  to  zero 
sample  mean  (spectra,  variances,  crest-to-trough  heights,  etc.).  Accord- 
ingly, there  was  no  particular  point  in  attempting  to  specially  consider 
biases,  and  it  was  decided  to  correct  ( t ) and  A^(t)  to  zero  sample  mean. 

The  programming  of  Eq.  (4)  and  (41)  is  straightforward.  The 
parameters  necessary  had  been  included  in  the  calibrated  data  files  des- 
cribed in  Reference  5.  The  total  radar  slant  range  was  approximated 
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according  to  the  assumptions  previously  mentioned,  and  two  new  time  histo- 
ries were  produced  for  each  recording  interval;  R^_  ( t ) and  Av(t).  One 

particular  interval  was  selected  for  the  purpose  of  detailed  numerical 
checking.  This  was  run  401  which  was  the  first  interval  on  analog  tape  145 
(index  17).  This  particular  run  was  chosen  because  the  waves  were  high, 
all  of  the  digitized  signals  were  relatively  large,  and  the  ship  was  in 
head  seas  at  low  speed  (8.7  knots). 

Figures  5a  and  5b  are  expanded  time  histories  covering  the  first 
300  odd  seconds  of  this  interval.  Only  the  5 data  channels  involved  in 
this  first  phase  of  the  data  reduction  are  shown.  At  the  top  is  the  slant 
(raw)  radar  range  corrected  to  zero  sample  mean  (dashed  line).  Superim- 
posed is  the  "corrected"  range  (Rc(t)),  also  corrected  to  zero  sample  mean. 
The  next  two  frames  down  each  page  are  scaled  roll  and  pitch.  Next  in  order 
are  the  "vertical"  accelerations.  The  solid  trace  is  A (t),  the  "true  verti- 
cal" acceleration.  Superimposed  in  dashed  lines  is  theVobserved  body 
vertical  acceleration,  A (t) . Finally  at  the  bottom  is  the  observed  body 
lateral  acceleration,  A^ft). 

In  the  process  of  checking  the  numerical  work  and  the  relationships 
between  the  various  channels  it  was  found  that  the  sense  of  pitch  as 
stated  in  Reference  5 was  opposite  to  what  would  be  expected.  At  the  rela- 
tively slow  ship  speed  involved  in  this  case  the  vertical  acceleration 
forward  should  be  nearer  in-phase  with  pitch  than  out  of  phase.  Also, 
the  derived  wave  elevations  (to  be  later  displayed)  implied  that  the  bow 
was  not  lifting  after  passage  of  a wave  crest.  The  strongest  evidence 
that  a sense  was  wrong  was  in  the  relationship  between  pitch  and  vertical 
acceleration.  Because  of  the  gravity  bias  in  the  vertical  acceleration 
signal  it  was  very  difficult  to  believe  that  the  sense  of  the  accelera- 
tion was  incorrect,  and  it  was  thus  concluded  that  the  pitch  sense  had  to 
be  reversed.  This  reversal  was  incorporated  in  the  programming  and  in 
Figure  5.  It  thus  should  also  be  noted  that  the  signs  of  the  pitch 
extremes  given  in  Reference  5 are  reversed. 

During  this  checking  process,  the  relationship  between  the  sample 
mean  of  pitch  and  electrical  zero  was  reviewed.  In  many  of  the  intervals 
digitized  from  the  second  season  the  sample  mean  differed  from  electrical 
zero  (the  stated  zero  pitch  condition)  by  a much  greater  amount  than  would 
be  expected  in  view  of  the  fact  that  (for  instance)  a 1/2°  pitch  bias 
corresponds  to  a 10  foot  ship  trim.  Accordingly,  it  was  felt  best  to 
correct  both  the  angular  responses  to  zero  sample  mean  before  processing. 

An  additional  observation  about  pitch  must  be  made.  It  may  be 
noted  that  the  pitch  trace  in  Figure  5 is  peculiarly  asymmetrical.  The 
large,  sudden  bow  down  (negative)  pitches  are  most  peculiar  in  view  of 
the  nominal  wave  conditions  and  ship  speed.  After  completing  the  deri- 
vation of  the  wave  trace  corresponding  to  Figure  5>  it  was  realized  that 
the  large  time  rates  of  change  of  bow  down  pitch  to  be  occasionally  seen 
apparently  corresponded  with  the  time  the  bow  was  plowing  into  particu- 
larly large  individual  waves.  It  is  thus  thought  possible  that  the 
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peculiar  asymmetric  negative  pitch  trace  results  from  transient  longi- 
tudinal deceleration  of  the  ship.  (A  sudden  longitudinal  deceleration 
has  the  effect  of  moving  the  pendulum  in  an  apparently  bow  down  direc- 
tion.) It  is  thus  probable  that  the  previously  made  zero  longitudinal 
acceleration  assumption  is  incorrect  at  times. 

The  object  of  this  phase  of  the  data  reduction  was  to  compute 
Rc(t)  and  ( t ) . It  is  obvious  from  Figure  5 that  the  difference  between 
A^(t)  and  A^(t)  is  exceedingly  small,  despite  the  fact  that  the  example 
involves  some  of  the  largest  motions  in  the  second  season  data  set.  Thus 
the  problems  about  pitch  just  noted,  and  the  previously  noted  question 
about  the  influence  of  sway  upon  the  roll  pendulum,  appear  not  to  be 
serious  with  respect  to  the  accelerat ions . 

Referring  again  to  Figure  5>  the  differences  between  slant  (raw) 
range  and  the  "corrected"  range,  Rc(t)  are  appreciable.  Analysis  of  the 
detail  of  the  numerical  result  revealed  that  the  dominant  terms  in  the 
correction  are  ^ , and  in  the  multiplier  of  Rg(t)  in  Equation 

(40).  In  Figure  5>  b3gQ,  is  effectively  the  roll  angle  times  half  the 
ship  beam.  Thus  for  the  cases  where  the  accelerometer  package  was 
located  in  the  radar  pedestal  the  corrections  should  be  much  less.  The 
reasons  for  the  dominance  of  the  terms  involving  b are  relatively  simple. 

3 2 q 

First,  the  actual  aiming  of  the  radar  relative  to  the  ship  was  at  90 
azimuth,  74  depression  --  almost  vertical.  This  in  conjunction  with  the 
cosine  nature  of  the  correction  makes  the  vertical  component  of  the  slant 
range  not  as  far  different  from  the  slant  range  as  it  could  have  been. 

(The  correction  is  mostly  just  the  cosine  of  the  roll  angle.)  Secondly, 
the  radar  package  in  the  case  of  Figure  5 is  located  nearly  abeam  of  the 
radar.  Thus  both  b^Q.^  amd  ( b^ ^ — 1 ) (Eq.40)  are  inherently  small.  The 
net  indication  from  this  is  that  the  problem  of  the  influence  of  sway  on 
the  roll  pendulum  is  still  a problem.  Though  the  problem  may  mostly 
influence  the  answer  for  the  voyages  where  the  accelerometers  were  not  in 
the  radar  pedestal,  it  was  thought  wise  to  plan  on  documenting  roll  in 
the  final  data  reduction  process. 


Correction  of  Radar  Range  for  Vertical  Motion  of  the  Ship 

The  purpose  of  the  work  is  to  produce  time  histories  of  wave 
el evat i on : 


-C(t)  = Rc(t)  + Zt(t) 


(42) 


and/or  the  encountered  spectrum  of  waves,  S^r(u)e).  Actually,  both  would 

be  desirable  results.  There  are  two  approaches  possible  to  make  an  effec- 
tive conversion  of  A^(t)  to  Z (t)  for  purposes  of  Eq.  (42). 


The  first  approach  involves  manipulations  of  the  acceleration  and 
range  spectra  and  cross-spectrum  to  result  in  a wave  spectrum.  In  this 
approach  rio  direct  estimate  of  -£(t)  is  made.  The  method  follows.  The 
spectrum  of  wave  elevation  (Eq.42)  may  be  written: 
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SCC(<"e)  = SRR(u)e)  + SZZ(tUe)  + 2CRZ(u)e) 


(43) 


where: 

SDD(u>  ) = spectrum  of  corrected  range,  R (t) 
kk  e c 

Szz(uie)  = spectrum  of  translation  2^(0 

CD7(o)  ) = co-spectrum  of  R (t)  and  Z.  (t) 

Kt  e c t 

Of  the  required  terms,  SDD(<u  ) may  be  estimated  since  R (t)  is  available 

KK  6 C 

from  the  previous  data  reduction  phase.  The  other  terms  cannot  be  esti- 
mated directly.  However  the  spectrum  of  A (t)  and  the  cross-spectrum  of 
R^_ ( t ) and  Av(t)  may  be  estimated: 

S..(tu  ) = spectrum  of  A (t) 

AA  e r v'  ' 

Cn.(u)  ) = co-spectrum  of  R (t)  and  A (t) 

RA  e r c v 

and  the  required  spectra  involving  Z(t)  may  then  be  estimated  as: 


S22(.e)  - 


SAA(*e> 


CRZ^ 


so  that: 


= SRR^e^  + x 4 a CRA^ 

e e 


So  long  as  the  estimated  spectra  derived  from  accelerations 
approach  or  can  be  made  to  approach  zero  sufficiently  fast  as  - 0 
there  are  no  particular  problems,  except  that  estimates  of  wave  eleva- 
tion in  the  time  domain  are  not  made. 


The  second  general  approach  is  to  perform  a numerical  double  inte- 
gration on  the  acceleration,  A^(t).  The  data  samples  are  of  20  minute 
duration  taken  every  four  hours.  Thus  no  matter  what  the  integration 
method  it  will  not  be  possible  to  produce  answers  for  fluctuations  which 
correspond  to  near  zero  frequency. 

Accordingly,  the  numerical  double  integration  may  be  thought  of 
as  band-pass  filtering,  just  as  in  some  forms  of  real  time  analog  inte- 
gration. Because  there  is  no  real-time  requirement  in  the  present  problem, 
the  double  integration  filter  response  function  does  not  have  to  be  real- 
izable, and  there  is  an  advantage  in  a non-real izable  filter  in  that 
incorrect  phase  lags  need  not  be  introduced. 

The  ideal  filter  transfer  function  for  double  integration  with 
respect  to  time  is:. 
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That  is;  the  phase  lag  is  180°,  independent  of  frequency,  and  the  ampli- 
tude response  is  ]/wJ . The  time  domain  equivalent  of  Eq.  (45)  is 
symmetrical  in  time;  that  is,  in  order  to  compute  the  instantaneous 
"present"  value  of  the  double  integral,  as  much  of  the  "future"  is 
required  as  of  the  "past."  This  is  not  a disadvantage  in  data  reduction 
since  everything  the  analyst  has  to  work  with  has  already  been  recorded  — 
at  any  instant  in  mid-sample  finite  amounts  of  "past"  and  "future"  are 
known. 


The  obvious  problem  in  Eq.  (45)  is  at  zero  frequency,  so  that  for 
any  sort  of  practical  implementation  the  transfer  function  specification 
must  be  altered  to: 


H<“e> 


r 


f(u.'e)  : 0 oue  < «JL'o 


-1/oj  ; x £ to 

e o e 


(46) 


For  reasons  previously  indicated  the  function  f(u>e)  can  be  zero  at  U)e=0 

to  make  the  filter  band-pass  and  must  grow  in  some  way  so  as  to  be 

— 1 / a)  ‘ at  in  . The  output  from  such  a band-pass  filter  will  be  correct 
o o 

double  integration  at  frequencies  exceeding  u)  and  quite  incorrect  inte- 
gration for  frequencies  below. 


The  Low  Frequency  Cut-Off 

In  the  initial  development  of  the  process  it  was  assumed  that  a 
fixed  low  frequency  cutoff  would  be  desirable.  The  accelerations  to  be 
processed  contain  the  ship's  rigid  body  motion.  The  only  moderately 
sharply  tuned  resonance  in  the  ship's  rigid  body  motion  is  that  of  roll. 
The  ship's  rolling  period  may  vary  between  20  and  30  seconds,  depending 
upon  loading.  It  must  be  assumed  that  even  after  correction  for  roll 
angle  there  will  be  at  least  some  accelerometer  output  content  in  this 
range  of  periods,  which  should  be  handled  correctly.  On  this  basis  the 
low  frequency  cutoff  should  correspond  to  encounter  periods  of  at  least 
60  seconds,  longer  if  possible. 

For  trial  purposes,  a value  of  0)q  of  0.0524  (120  second  encounter 
period)  was  assumed  and  the  implications  with  respect  to  double  integra- 
tion errors  were  studied  in  the  wave  frequency/heading  plane.  Essential 
results  are  contained  in  Figures  6 and  7.  These  are  polar  plots,  wave 
frequency  is  along  the  radius  and  an  auxiliary  wave  length  scale  is  shown. 
The  regions  are  shown  for  combinations  of  wave  length  and  heading  which 
produce  roll  resonance  for  roll  periods  between  20  and  30  seconds. 

A second  pair  of  regions  is  indicated  where  encounter  periods  are  greater 
than  120  sec.  These  last  are  the  regions  where  double  integration  errors 
will  exist  for  uu  = .0524,  Eq.  (46). 
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The  small  region  of  integration  error  around  the  origin  is  about 
the  same  size  for  either  20  or  30  knot  ship  speeds.  The  wave  lengths 
involved  are  measured  in  tens  of  thousands  of  feet.  According  to  con- 
ventional thinking,  wave  components  of  this  length  either  do  not  exist 
or  are  of  such  small  amplitude  that  they  cannot  be  resolved,  and  thus 
the  ship  motions  response  to  such  wavelengths  may  be  assumed  negligible. 

On  this  basis  no  significant  influence  on  the  double  integration  of  the 
low  frequency  cutoff  may  be  anticipated  in  head  through  beam  seas  to  a 
wave  heading  of  around  70°. 

The  second  region  of  integration  error  is  centered  around  the  line 

(shown  dashed)  corresponding  to  quartering  or  following  wave  components 

which  travel  at  the  ship  speed.  Along  the  dashed  line  in  Figures  6 and  7 

the  encounter  period  is  infinite  (ui^=0) . A proper  double  integration  is 

not  possible  with  discontinuous  sample  records.  The  position  of  the 

dashed  line  is  fixed  by  ship  speed.  The  width  of  the  error  region  is 

fixed  by  ujq.  Considering  the  available  experimental  and  theoretical  ship 

motions  results,  significant  ship  motions  (pitch,  heave,  roll)  are  not 

expected  for  waves  shorter  than  about  1/4  ship  length  (say  200  feet  in  the 

present  case)  unless  roll  resonance  is  involved.  Considering  Figure  6, 

the  second  error  region  corresponds  to  waves  shorter  than  200  feet  for 

quartering  headings  to  25°,  and  at  worst  to  waves  shorter  than  250  feet 

for  following  seas.  Serious  integration  errors  are  thus  not  expected  for 

this  ship  speed.  Similarly,  Figure  7 indicates  that  integration  errors 

should  not  be  significant  over  headings  from  head  to  almost  quartering 

seas  (50°)  for  ship  speeds  to  30  knots.  For  following/quartering  seas 

(0  to  50°  heading)  integration  errors  will  be  introduced  for  wave  lengths 

which  should  produce  ship  motion.  This  type  of  error  cannot  be  eliminated, 

its  magnitude  can  only  be  reduced  by  decreasing  u)  . 

o 


Figure  8 is  an  enlargement  of  the  quartering/following  sea  error 
regions.  In  this  case  the  radius  of  the  plot  is  wave  length.  In  the 
20  to  30  knot  operating  speed  range  of  the  present  ship,  errors  can  be 
expected  which  correspond  to  the  very  significant  (for  the  Atlantic) 
range  of  300  to  600  feet  wave  lengths. 

It  seemed  reasonably  clear  from  the  foregoing  that  the  low-frequency 
cutoff,  u)q  (Eq.46)  should  not  exceed  0.0524  rad/sec  in  the  present  applica- 
tion. The  widths  of  the  bands  in  Figure  8 are  roughly  proportional  to  u i 
so  that  the  influence  of  the  integration  errors  for  following  seas  should 
reduce  in  proportion.  However  there  seemed  to  be  little  point  in  exerting 
large  effort  to  reduce  the  cutoff  frequency  much  beyond  1/2  of  that  assumed 
for  Figures  6-8.  For  wave  components  of  near  zero  encounter  frequency  the 
slant  range  will  also  vary  slowly,  so  that  the  apparent  period  of  the  wave 
components  will  approach  the  sample  length.  In  this  circumstance  there 
will  not  be  anough  information  available  to  make  reasonable  statistical 
estimates  about  the  component  wave  lengths  in  question.  Nothing  possible 
in  analysis  will  help. 


31 


Considering  the  sample  length  of  20  minutes,  component  wave 
encounter  periods  of  120  to  400  seconds  (ui^  = 0.0157  to  .0524)  involve 
between  3 and  10  encounter  cycles  in  the  sample.  It  therefore  seemed  a 
reasonable  compromise  for  numerical  development  purposes  to  set  a target 
cut-off  of  about  0.0262  rad/sec  (encounter  period  240  sec). 


Filtering  Methods  - Sampled  Data 


By  taking  the  Fourier  transform  of  Eq.  (46)  a series  of  weighting 
coefficients  may  be  determined  so  that  filtering  may  be  done  with  a dis- 
crete convolution  of  the  form: 


Z(rAt) 

where 


p 

C A (nAt)  + I C.[a  (nAt-jAt)  + A (nAt+jAt)] 

ov  ..IV  J v J 

j = l 

At  = sampling  interval  of  time  series 


C. 

J 


weighting  coefficients  for 
non-real izable  filter. 


(47) 


In  a sample  of  N points,  n=l,2...N,  there  is  a starting  and  stop- 
ping transient  (P)  points  in  length  which  is  only  partially  correct. 

(To  compute  Z(0),  (P)  data  points  for  A (nAt)  are  required  for  negative 

n --  these  points  are  not  known.)  Accordingly,  it  is  prudent  to  compute 
Z(nAt)  only  for  (P+1)  s n s (N-P) . This  has  two  results: 

1.  Results  are  obtained  for  a shorter  sample.  In  effect 
the  resulting  time  series  for  Z(t)  is  2P  points  shorter 
than  that  of  the  original  (A^(t)). 

2.  The  computing  time  is  approximately  proportional  to  the 
njmber  of  multiplications  which  are  approximately  P(N-2P). 

An  alternate  approach  is  the  use  of  "fast  convolution"  methods 
using  the  FFT  algorithm.  Broadly,  in  this  approach  the  direct  discrete 
Fourier  transform  of  the  data  is  taken.  This  results  in  a complex 
"spectrum"  of  the  original  data.  The  filtering  is  done  in  the  frequency 
domain  by  multiplying  this  "spectrum"  by  the  fi'ter  transfer  function. 
After  this  operation  the  inverse  discrete  Fourier  transfor  is  taken  and 
the  result  is  the  filtered  time  series  in  the  time  domain.  The  discrete 
Fourier  transform  has  as  an  underlying  assumption  that  the  data  Is  peri- 
odic in  the  sample  length.  In  the  present  case  this  is  not  true  so  that 
an  additional  pre-processing  step  is  necessary  to  mini  ize  end  effects. 
This  takes  the  form  of  "tapering"  or  attenuating  the  ends  of  the  data  so 
that  the  first  and  last  points  are  near  zero.  After  the  direct  trans- 
form, multiplication,  and  inverse  transforming,  the  result  is  a "tapered" 
rendition  of  the  filtered  data.  Accordingly  data  at  the  ends  of  the 
record  are  also  lost  in  the  fast  convolution  method. 

For  comparison  of  methods  the  FFT  method  may  be  hypothesized  in 
the  following  steps. 
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1.  Pre-process  data  with  the  Tukey  10%  cosine  taper.  This  is 
a standard  tapering  used  in  FFT  spectrum  analysis.  10%  of 
the  data  at  each  end  of  the  sample  is  tapered  so  that  the 
final  results  will  be  a time  series  N/5  points  shorter  than 
the  original  N point  series.  Approximate  number  of  multi- 
pl i cat i ons  = N/5 . 

2.  Do  Direct  FFT.  Approximate  number  of  real  multiplications  = 2MN 
where  M is  defined  by  N=2  . 

3.  Multiply  direct  transform  by  transfer  function  --  2N  multi- 
pl i cat i ons . 

4.  Do  Inverse  FFT  — 2MN  multiplications. 

Total  number  of  multiplications  for  this  process  is  (2.2+4M)N. 

The  foregoing  provided  some  basis  on  which  to  make  a choice  between 
the  methods.  From  previous  experience  with  digital  band-pass  filter  coef- 
ficients of  the  form  of  Eq.  (47)  it  was  estimated  that  the  minimum  value 
of  (P)  for  a reasonably  precise  filtering  job  would  be  about  equal  to  the 
number  of  points  necessary  to  define  in  the  time  domain,  the  period  corres- 
ponding to  the  peak  of  the  transfer  function.  Assuming  the  minimum  accept- 
able period  to  be  120  sec  or  1/10  of  sample  length  (N),  it  appeared  that 
(P)  should  be  about  800,  since  the  sample  length  for  the  present  data  is 
8192. 

It  thus  appeared  that  the  data  loss  would  be  the  same  in  either 
method  — approximately  20%  of  the  sample  length  — and  no  choice  between 
the  methods  was  dictated. 

Considering  computer  time  to  be  proportional  to  the  number  of  real 
multiplications,  the  estimated  ratio  of  computation  times  may  be  formed: 

w _ Computation  Time  for  Eg,  (47) 

Computation  Tine  for  Fast  Convolution  Method 

P(N-2P) 

545  (2. 2t4m)  N 

Noting  that  P 800  (minimum)  and  M= 1 3 in  the  present  program: 

W > 12 

Thus  the  fast  convolution  method  appeared  to  be  an  order  of  magnitude 
faster  in  the  present  application.  Because  of  the  large  amount  of  data 
involved  this  speed  advantage  was  of  practical  significance. 

The  fast  convolution  method  was  clearly  the  best  choice.  It  had 
the  further  advantage  that  mid-way  through  the  process  nearly  half  of 
the  work  involved  in  obtaining  the  spectrum  and  cross  spectrum  of  the 
vertical  translation  and  corrected  slant  range  (Eq.43)  would  be  done  and 
could  be  salvaged. 
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Preliminary  Numerical  Development;  Fast  Convolution  Method 

Initial  trials  with  the  fast  convolution  method  were  carried  out 
with  cosine  wave  inputs.  Essentially  a cosine  wave  of  frequency  uj  was 
assumed  to  be  the  acceleration  and  a time  series  of  cos (gj  jAt)  was  gen- 

/ \ O 2 

erated  for  j=0, 1 . . . (N- 1) . This  series  was  multiplied  by  ~'/^0  to 
generate  a "true  displacement"  series.  The  processing  was  carried  out 
on  the  "acceleration"  series  and  the  results  were  compared  with  the  "true 
displacement"  series  to  assess  precision. 

If  the  sample  length  is  an  exact  multiple  of  the  periodicity  of 
the  data,  the  simple  process  of  taking  the  inverse  Fourier  transform  of 
(-1/a)  times  the  Direct  Fourier  Transform)  was  found  to  be  nearly  as 
precise  as  the  number  of  significant  digits  carried  in  the  computation. 
However  the  data  in  question  has  no  such  restrictions  about  periodicity. 
The  simple  process  just  outlined  did  not  result  in  spectacularly  good 
precision  for  input  periods  not  commensurate  with  sample  length.  In 
this  case  the  fractional  cycle  at  the  end  of  the  sample  duration  causes 
the  Fourier  Spectrum  to  be  higher  than  it  should  be  throughout  the  fre- 
quency range.  Multiplication  of  the  direct  transform  by  -\/x  thus 
unduly  emphasizes  any  spurious  low  frequency  Fourier  spectrum  components 
because  the  spurious  components  arise  through  "leakage"  of  variance  from 
the  actual  pure  input. 

With  real  data  the  exact  frequency  content  is  not  known  in  advance 
and  it  should  not  be  assumed  that  the  data  is  free  of  linear  trends  since 
an  appreciable  linear  trend  also  has  the  effect  of  introducing  spurious 
leakage  from  high  to  low  frequencies.  Accordingly  the  simple  process 
must  be  augmented  by  two  pre-processing  steps; 

1.  Correct  sample  to  zero  mean  and  remove  linear  trends  in 
one  of  the  standard  manners.  In  the  present  case  the 
method  chosen  was  the  minimum  computation  method  recom- 
mended by  Blackman  and  Tukey  for  spectrum  analysis. 

2.  Perform  the  1 0%  cosine  taper  as  recommended  by  Tukey  to 
reduce  "leakage."  This  procedure  results  in  an  effective 
discard  of  10%  of  the  final  data  at  each  end  of  the  series 
as  mentioned  in  a previous  section.  It  also  is  one  of  the 
required  steps  for  FFT  spectrum  analysis. 

The  i ncorporat ion  of  these  two  steps  resulted  in  about  the  same 
precision  of  displacement  estimates  whether  or  not  the  acceleration  input 
period  was  commensurate  with  sample  duration.  (For  commensurate  input 
periods  the  precision  was  degraded,  for  incommensurate  periods  it  was 
improved.)  The  magnitude  of  rms  prediction  error  was  about  2 or  3%  of 
displacement  amplitude  for  periods  the  equivalent  of  60  seconds  real  time. 
However  the  errors  did  not  fall  off  as  rapidly  as  was  expected  as  period 
was  shortened.  With  some  simple  trial  and  error  experimentation  the  pre- 
cision was  found  to  be  a function  of  the  cutoff  uj  , Eq.  (46),  and,  to  a 
much  smaller  degree,  of  the  shape  of  f (u>e) , Eq . (46). 
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The  reasons  are  purely  numerical.  If  the  input  acceleration  is  of 
very  high  frequency  the  FFT  acceleration  spectrum  has  essentially  noise  at 
low  frequencies.  If  the  input  frequency  is  high  enough  the  1/u)  charac- 
teristic of  the  filter  can  result  in  the  spurious  low  frequency  components 
of  the  acceleration  spectrum  being  transformed  into  relatively  significant 
components  of  the  displacement  spectrum.  In  addition,  with  tapering,  the 
spectrum  estimates  at  frequencies  just  below  those  corresponding  to  the 
lowest  frequency  actually  present  (oj.,  say)  correspond  mostly  to  leakage 

from  the  lowest  frequency  component.  The  appropriate  double  integration 
multiplier  for  these  near  frequencies  may  consequently  be  (-l/u).4 S)  rather 
than  the  reciprocal  of  the  square  of  the  actual  frequency  involved. 

These  preliminary  results  with  single  frequency  accelerations 
suggested  that  a fixed  low  frequency  cut-off  characteristic,  as  was  orig- 
inally envisioned,  would  not  be  the  optimum  solution  for  the  present  work. 
Roughly  half  of  the  present  data  is  in  nominal  head  or  bow  seas,  and  the 
other  half  is  in  nominal  following  or  quartering  seas.  One  constant  low 
frequency  cutoff  suitable  for  both  situations  did  not  seem  possible. 


Radar  Wave  Processing  System 

It  was  clear  that  in  order  to  make  a best  estimate  of  doubly  inte- 
grated accel erat ion,  the  acceleration  spectra  would  have  to  be  developed 
simultaneously  in  order  to  make  a best  choice  of  low  frequency  cutoff. 

Thus  what  initially  appeared  attractive  to  minimize  computation  time 
became  a requirement.  Accordingly  the  required  (largely  standard)  sub- 
programming was  gathered,  the  system  was  checked  out  with  simulated  data, 
and  by  this  means  a final  form  was  developed  for  the  double  integration 
filter.  All  of  the  individual  standard  steps  in  the  process  are  described 
in  Reference  10. 

The  steps  in  the  computation  are  as  follows  for  each  interval  which 
was  processed.  (All  reading  and  storing  operations  involved  disk  files 
accessible  by  the  computer.) 

1.  Read  the  basic  digitized  data  (slant  range,  roll,  pitch, 
vertical  and  lateral  accelerations)  calibrate  each  channel, 
and  compute  R^(t)  and  A^(t)  in  accordance  with  the  previous 
discussion.  These  new  time  series  were  placed  in  8192  point 
arrays,  the  range  in  array  XR  and  the  acceleration  in  array  XI. 

2.  Correct  the  two  time  series  ( X R , X I ) for  linear  trends  and  to 
zero  mean  in  accordance  with  the  minimum  computation  method 
recommended  by  Blackman  and  Tukey  (Ref. 10). 

3.  Store  the  range,  XR  in  a temporary  file  for  later  use.  (When 
the  double  integration  is  completed  the  range  must  be  retrieved 
in  order  to  compute  wave  elevation.) 

4.  Apply  the  10%  Tukey  cosine  taper  (Ref. 10)  to  both  the  range  and 

acceleration  arrays. 
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5.  Compute  the  direct  fast  Fourier  Transform  (Ref. 10)  of  arrays  XR 
and  XI,  assuming  that  array  XR  (range)  is  the  real  part  and 
array  XI  (acceleration)  is  the  imaginary  part  of  a complex 
function  of  time.  (This  approach  allows  spectra  and  cross 
spectra  to  be  developed  later  with  about  a 60%  saving  in  com- 
putation time  over  an  approach  in  which  the  spectra  were 
developed  in  a separate  operation.)  The  results  of  this  step 
is  a scrambled  transform  in  arrays  XR  and  XI  --  unscrambling 

is  a matter  of  simple  arithmetic. 

6.  Store  both  parts  of  the  direct  transform  in  a second  temp- 
orary file  for  later  use. 

7.  Compute  over-resolved  scalar  spectra  of  range  and  accelera- 
tion from  the  transform,  placing  the  results  in  arrays  XR 
and  XI  in  the  process.  These  spectra  are  defined  at  the 
minimum  delta  frequency  possible  with  the  given  size  arrays 
(Au)  = 2^/8192  At  Rs  .005  rad/sec). 

8.  From  the  over-resolved  acceleration  spectrum  compute  the 
parameters  (to  be  described)  of  the  double  integration 

f i Iter. 

9.  Read  back  the  scrambled  direct  transform  ( XR , X I ) stored  in 
step  6. 

10.  Multiply  the  acceleration  transform  by  the  double  integra- 
tion filter.  This  was  accomplished  by  unscrambling  the 
range  and  acceleration  transforms,  multiplying  the  accel- 
eration transform  by  the  filter  response,  and  then 
re-combining  the  two  transforms  so  that  the  result  (in 
arrays  XR  and  XI)  was  as  though  the  direct  Fast  Fourier 
Transform  had  been  done  with  arrays  containing  range  and 

d i s p 1 acement  rather  than  range  and  acceleration. 

11.  From  the  results  of  step  10  compute  the  scalar  spectra  of 
range  and  displacement  and  the  range-d isplacement  cross- 
spectrum. In  this  step  arrays  XR  and  XI  were  not  modified. 

The  spectra  and  cross-spectrum  could  be  placed  in  other 
smaller  arrays  at  this  stage  because  the  final  averaging 
over  blocks  of  basic  estimates  (Ref. 10)  was  carried  out. 

12.  Using  the  results  of  step  11,  develop  the  radar  wave 
spectrum  in  accordance  with  Eq.  (43),  retaining  for  record 
purposes  the  rms  corrected  range  R (t)  and  the  rms  dis- 
placement . 

13.  Replace  arrays  XR  and  XI  by  the  transform  of  displacement 
alone  in  an  unscrambling  step  similar  to  that  performed 
previously.  The  result  is  as  though  the  direct  FFT  had 
been  done  with  displacement  in  array  XR  and  array  XI 
zeroed . 

14.  Compute  the  inverse  Fast  Fourier  Transform  of  (XR+iXI). 

The  result  is  the  (tapered)  time  history  of  vertical 
displacement  in  array  XR. 
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15.  Store  the  displacement  time  history  for  later  use. 

(initially  this  was  a de-bugging  step,  but  it  was  retained 
later  for  use  in  processing  the  Tucker  Meter  data.) 

16.  Read  the  range  R (t)  from  the  first  temporary  file  (step  3) 
into  array  XI . 

17.  Compute  the  wave  elevation  in  accordance  with  Eq.  (42)  as: 

C(t)  = -Rfl(t)  - Zt(t) 

= -XR  - XI 

18.  Write  out  a summary  file  of  results  for  the  interval. 

Simulation  of  Data  for  Final  Development  of  Low  Frequency 

Cutoff  Parameters  of  H(ou  ) 

e 

In  order  to  develop  a reasonable  system  for  adapting  the  low  fre- 
quency cutoff  parameters  of  H(u)g)  to  the  acceleration  data,  a more 
realistic  simulation  of  input  data  was  in  order.  Because  a simulation  of 
both  the  acceleration  and  the  corresponding  displacement  was  required,  the 
approach  taken  was  to  sum  a number  of  co-sinusoids  with  incommensurate 
frequencies,  equal  amplitudes  and  random  starting  phases.  Accordingly, 
the  "true"  displacement  may  be  computed  from  the  acceleration  components. 
The  simulated  acceleration  was  processed.  The  simulated  displacement  was 
used  for  comparison  with  the  results  to  assess  precision. 

In  particular  the  simulation  was  as  follows: 

NS 

A^(t)  = a S cos(ujjt  + Gj)  (48) 

NS 

Z (t)  = -a  I (lAu.s)cos(uj.t  + e.)  (49) 

1 j=l  J J J 

where  a = /2/NS 

= desired  rms  acceleration 
NS  = number  of  sinusoids  summed 

e.  = random  initial  phase 

(uniform  in  0 s s.  s 2n) 

J 

u).  = frequencies. 

J 

To  derive  the  frequencies,  u).,  an  acceleration  spectrum  was  assumed 
in  the  form: 
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= 0 elsewhere 


The  spectrum  was  broken  into  (NS+1)  equi-energy  bands  and  the  frequencies 
iu.  were  determined  by  solving: 

0). 

. 2 -* 

J'SAWd“ 

for  j=l ,2. . .NS 

This  procedure  substantially  insures  incommensurate  frequencies.  The 
result,  in  Eq.  (48),  for  NS=1  is  a pure  co-sinusoid  of  frequency  (tu1+u^)/2. 
As  NS  is  increased  the  result  begins  to  resemble  a band-passed  random 
funct ion. 

In  virtually  all  the  simulations  using  the  system  outlined  in  the 
last  section  white-noise  slant  range  "data"  with  a standard  deviation  of 
10  times  that  of  the  simulated  acceleration  was  also  simulated.  Some 
trials  were  also  made  with  zeroed  simulated  slant  range.  There  appeared 
to  be  no  numerical  cross-talk  between  the  slant  range  and  acceleration. 

Final  Version  of  the  Low  Frequency  Cutoff  Parameters 

The  final  version  of  the  low-frequency  cutoff  parameters  was 
developed  more  or  less  by  trial  and  error.  Harmonic  acceleration  and  both 
narrow  and  wide-band  "random"  accelerations  were  processed.  The  results 
were  compared  with  the  simulated  true  displacement  and  the  precision  was 
assessed  by  computing  the  rms  error  over  the  middle  75%  of  the  sample  to 
eliminate  the  distortions  induced  by  "tapering."  This  rms  error  was 
divided  by  the  rms  simulated  displacement  and  the  result,  in  percent, 
formed  the  basis  for  judgment  of  the  adequacy  of  the  process. 

Figure  9 summarizes  the  final  computation  system  chosen  for  H(uj  ). 

8 6 

Given  the  acceleration  spectrum,  the  acceleration  variance,  a ^ , is  known. 

A measure  of  the  lower  limit  of  significant  acceleration  spectral  density 
is  formed  by  finding  a frequency,  id  such  that  2%  of  total  variance 

is  attributable  to  lower  frequencies.  This  approach  was  chosen  to  avoid 
odd,  possibly  spurious,  spikes  in  an  over-resolved  analysis  from  upset- 
ting things. 

The  cutoff  frequency,  Eq.  (46),  is  formed  by  taking  75 % of  u:  qm 

and  adding  an  empirically  determined  term  (Fig. 9)  which  has  a negligible 
effect  on  the  cutoff  for  values  of  id  which  are  high  relative  to  the 

folding  frequency.  The  term  significantly  alters  the  cutoff  only  for 
frequencies  corresponding  to  fewer  than  20  periods  in  the  sample  length. 
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6Ja  = 0.75"  CQ oa  -+-  3-5*2  /» 

^=2  cOo/3 

FIGURE  9 DEFINITION  OF  H(u)e) 
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The  shape  of  the  resulting  H(u>  ) is  indicated  in  Figure  9.  The 
2 ® 

function  is  ( — 1 / a>  ) above  the  cutoff,  m . It  is  constant  and  equal  to 
2 ^ 

(-1/uuq  ) from  (2u)q/3)  to  u)Q,  and  decreases  linearly  to  zero  for  lower 

frequencies.  It  is  thus  a rough  approximation  to  a band-pass  filter  as 
required.  Since  all  functions  are  defined  in  discrete  steps  in  frequency, 
the  actual  values  of  cutoff  frequencies  used  are  slightly  different  than 
shown  in  Figure  9 by  the  effects  of  truncation  in  transforming  frequencies 
to  the  corresponding  integer  steps. 

Results  of  Final  Precision  Checks 

The  results  of  the  processing  of  simulated  harmonic  and  narrow  band 
accelerations  are  shown  in  Figure  10.  The  frequency  abscissa  includes  the 
entire  range  of  the  potential  data.  The  chosen  folding  frequency  is 
3.33  Hz,  the  probable  maximum  resolvable  wave  encounter  frequency  is  about 
1.4  Hz.  The  narrow  band  excitation  is  indicated  by  circles  with  a hori- 
zontal line.  The  horizontal  line  indicates  the  bandwidth  (0^,1^  in  the 

section  describing  the  simulation).  Within  this  band  20  cosinusoids  were 
summed  in  random  phase  to  produce  the  acceleration. 

The  ordinate  is  rms  prediction  error  in  percent  of  rms  true  dis- 
placement. Simulations  were  carried  out  for  frequencies  from  0.08  to 
0.7  Hz  also.  The  rms  error  in  these  cases  was  less  than  0.01%,  and  was 
probably  controlled  by  the  number  of  significant  digits  (7)  carried  in 
the  computation. 

In  interpreting  the  results  and  developing  the  cutoff  parameters, 
errors  less  than  1%  were  considered  good,  and  the  basic  problem  was  to 
reduce  prediction  error  to  the  order  of  2%  or  less  for  frequencies  of 
0.008  Hz  (120  sec  encounter  periods).  Figure  10  represents  about  the 
best  that  could  be  done  in  this  respect.  It  became  clear  that  the  previ- 
ously mentioned  target  of  a system  good  down  to  a 240  second  encounter 
period  could  not  be  met.  Figure  10  indicates  that  the  basic  procedure  is 
good  for  encounter  periods  as  long  as  60  seconds  and  probably  acceptable 
from  60  to  120  second  encounter  periods.  For  longer  encounter  periods 
larger  errors  develop.  Their  importance  is  outlined  in  Figures  6-8  and 
a previous  section. 

Figure  11  contains  the  results  of  the  simulation  for  wide-band 
acceleration.  The  width  of  the  acceleration  spectrum  is  indicated  by 
the  sketch  of  a band-passed  spectrum  which  serves  as  a plotting  symbol. 

Each  wide-band  simulation  was  produced  by  summing  200  cosinusoids.  Rela- 
tive to  Figure  10,  the  ordinate  of  Figure  11  is  shifted  by  one  decade  so 
that  the  rms  error  for  all  cases  is  1%  or  less  of  true  rms  displacement. 

The  trend  line  from  Figure  10  is  shown.  The  magnitude  of  rms  error  is 
determined  largely  by  the  frequency  of  lowest  significant  spectral  density. 

In  all  the  previous  development  the  acceleration  data  input  to  the 
system  had  seven  digit  significance.  This  is  roughly  analogous  to  an 
analog  instrument  with  a dynamic  range  greater  than  150  db.  In  the  present 
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case  there  is  nothing  like  this  dynamic  range  inherent  anywhere  in  the 
system  until  the  final  data  processing  stage. 


The  effect  of  reduced  significance  was  studied  by  artificially 
reducing  the  resolution  of  the  simulated  acceleration.  For  instance, 
to  reduce  the  resolution  to  p%  of  rms  acceleration  a 

1.  The  simulated  acceleration  was  divided  by  .Olpc^. 

2.  This  result  was  truncated  by  discarding  digits  to  the 
right  of  the  decimal  point. 

3.  Finally  the  last  result  was  multiplied  by  O.Olpc^  to 
restore  scale. 

Many  of  the  simulations  shown  in  Figures  10  and  II  were  repeated 
with  input  data  of  various  levels  of  significance.  The  results  were 
remarkably  consistent  regardless  of  the  acceleration  bandwidth  or  fre- 
quency. The  estimated  precision  of  the  double  integration  process 
involving  real  data  may  be  summarized  as  follows: 

For  input  acceleration  data  having  an  rms  of  and  which  is 
resolved  to  p%  of  the  precision  of  the  prediction  of  time  domain 
displacement  in  the  middle  75%  of  the  sample  appears  to  be  as  indicated 
in  Figures  10  and  II,  £r  equal  to  (.6  to  1.0) (p),  whichever  is  greater 
for  the  frequency  content  involved. 

In  other  words,  if  the  acceleration  resolution  is  2%  of  rms 
acceleration,  the  rms  error  in  displacement  prediction  is  expected  to 
be  slightly  less  than  2?/  of  rms  displacement  unless  the  frequencies 
involved  are  less  than  0.01  Hz. 

In  digitizing  analog  magnetic  tape  the  ±1.5  volt  full  scale  output 
is  (for  practical  reasons  associated  with  existing  hardware)  resolved  to 
about  1 part  in  620  or  about  0.16%  of  full  scale.  For  most  random  pro- 
cesses, if  the  rms  is  in  excess  of  1/3  full  scale,  substantial  saturation 
may  be  involved.  Accordingly  under  the  best  possible  conditions  the  pre- 
cision of  the  time  domain  estimates  of  displacement  will  not  be  better 
than  about  1/2%. 


Sample  Evaluation  with  Actual  Data 

As  in  the  development  of  the  angular  corrections  to  range  and 
acceleration,  extensive  checking  of  the  process  was  carried  out  for  one 
interval.  The  interval  was  the  same  as  that  previously  chosen,  Run  401. 
Figures  12a  and  12b  are  sample  time  histories  obtained  from  the  reduction 
of  this  interval.  The  traces  at  the  top  of  the  figure  are  the  corrected 
range,  Rc(t),  and  the  corrected  acceleration,  Av(t).  Both  traces  have 
previously  been  presented  in  Figure  5.  The  "vertical  displacement"  in 
Figure  12  is  the  double  integration  of  the  acceleration,  and  the  trace 
marked  "wave"  is  the  final  time  domain  answer,  £(t). 
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The  displacement  trace  in  Figure  12a  displays  the  peculiarity 
brought  on  by  tapering.  The  samples  are  20-1/2  minutes  long  and  the 
cosine  taper  is  applied  to  the  first  and  last  10%  of  the  sample.  Thus 
approximately  the  first  and  last  120  seconds  of  the  displacement  result 
are  intentionally  distorted.  Only  the  last  part  of  the  displacement  in 
Figure  12a  is  considered  valid.  Thus  the  wave  trace  (derived  in  part 
from  the  displacement)  must  be  considered  invalid  for  the  first  and  last 
2 minutes  of  the  record  so  that  the  wave  time  histories  resulting  from 
the  data  reduction  process  are  at  best  valid  only  in  the  middle  16-1/2 
minutes  of  the  sample. 


The  wave  trace  displays  the  noise  which  is  inherent  in  the  cor- 
rected range.  Most  of  this  noise  seems  to  have  been  in  the  original 
slant  range  (Fig. 5),  though  some  must  certainly  have  been  induced  by 
the  correction  process.  A.  few  of  the  secondary  bumps  in  the  wave  trace 
border  on  the  hard-to-bel i eve.  While  it  is  felt  that  they  are  physically 
possible  for  the  severe  conditions  involved  in  the  example  it  must  be 
admitted  that  they  could  be  the  result  of  a combination  of  the  types  of 
errors  previously  described. 

Figure  13  has  been  prepared  to  illustrate  the  frequency  domain 
side  of  the  data  reduction  procedure  for  the  trial  interval,  Run  401. 

At  the  top  of  the  figure  are  plotted  the  spectra  of  corrected  range, 

R (t),  and  of  the  doubly  integrated  acceleration,  Z (t) . Also  indicated 

3 t 

is  the  cutoff  frequency,  u)q,  for  the  double  integration  filter.  In  the 
example  virtually  all  of  the  range  variance  is  associated  with  frequencies 
above  the  cutoff.  From  these  spectra  the  significant  double  amplitude  of 
radar  range  (4  x rms)  is  50.3  feet,  that  for  the  vertical  displacement  of 
the  radar  is  42.8  feet.  The  radar  is  moving  up  and  down  almost  as  much  as 
the  change  in  its  indicated  range.  In  the  middle  of  the  figure  is  shown 
the  sum  of  the  displacement  spectrum  and  twice  the  range-displacement 
co-spectrum.  This  is  effectively  an  additive  correction  to  the  range 
spectrum  in  accordance  with  Eq.  (43).  At  the  bottom  of  the  figure  is  the 
resulting  wave  spectrum.  The  correction  is  quite  large.  The  significant 
double  amplitude  of  wave  is  estimated  from  the  spectrum  to  be  36.0  feet. 
For  frequencies  above  0.45  the  correction  indicates  that  range  and  dis- 
placement are  largely  out  of  phase  so  that  the  net  (the  wave)  is  much 
smaller  than  the  corrected  radar  range.  Below  a frequency  of  0.45  the 
correction  is  positive,  which  means  that  the  radar  was  moving  up  and  down 
with  the  local  wave  so  that  the  radar  range  is  an  underestimate.  If  it 
is  assumed  that  the  ship  is  in  long-crested  head  seas  the  frequency  range 
where  the  correction  is  positive  corresponds  to  waves  roughly  1-1/4  to 
1-1/2  the  ship  length,  and  the  indication  is  plausible  that  radar  and 
local  wave  are  somewhat  in-phase. 

After  all  the  detailed  numerical  checking  and  the  analysis  of 
sample  results  for  internal  plausibility,  there  is  always  the  question 
about  whether  or  not  the  final  answer  is  in  the  "ball  park."  In  the  case 
of  the  sample  interval  the  log-book  data  gives  an  estimated  wave  height 
of  25  feet  and  a swell  height  of  20  feet.  The  raw  Tucker  Meter  wave  esti- 
mate was  a 13  foot  significant  height  (4  rms).  These  should  be  compared 
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with  the  radar  wave  estimate  of  36  feet  — that  is,  the  radar  indication 
was  nearly  three  times  that  of  the  Tucker  meter,  and  something  like 
1-1/2  times  the  visual  estimates. 

On  the  subject  of  visual  estimates,  the  interval  in  question  was 
recorded  at  2000  GMT,  10  January  1974,  at  a position  roughly  43°N  by 
24°  West,  in  what  was  described  as  overcast  weather  with  rain.  Local 
time  would  have  been  1800  hours  — probably  dark  or  at  least  twilight. 

The  wave  height  at  1000  hours  the  same  day  was  described  as  30  feet  with 
15  foot  swell,  and  that  the  following  morning  as  30-35  feet  with  20  foot 
swell.  During  this  period  of  nearly  a day  the  ship  speed  was  fairly  con- 
sistently 12  knots  or  less  except  for  an  apparently  brief  period  of  trying 
an  increased  speed.  It  accordingly  seemed  permissible  to  speculate  that 
the  visual  estimates  for  the  interval  in  question  could  have  been  on  the 
low  side.  If  the  waves  during  the  night  were  not  much  different  than  in 
the  afternoon  preceding  and  the  morning  following,  the  comparison  between 
visual  and  radar  would  have  been  30  as  opposed  to  36  feet,  not  too  bad 
a ballpark  estimate. 

In  advance  of  the  analysis  of  a large  number  of  runs  the  above 
speculations  were  not  altogether  reassuring.  One  other  simple  ballpark 
comparison  involves  computing  the  midship  longitudinal  stress  spectrum, 
dividing  this  by  the  wave  spectrum  derived  from  the  radar,  and  taking  the 
square  root.  The  result  of  this  operation  (where  both  spectra  are  appre- 
ciable) is  an  estimate  of  the  effective  stress  response  amplitude  per  foot 
of  wave  amplitude.  The  stress  spectrum  for  Run  401  was  computed  and  the 
indicated  operation  carried  out.  A nominal  midship  section  modulus  was 
found  from  Reference  3 and  this  used  to  convert  the  result  to  a moment 
response  in  (foot  tons/foot).  The  results  are  shown  in  Figure  14.  The 
abscissa  on  the  plot  is  wave  frequency.  Wave  frequency  plotting  positions 
for  the  estimates  from  the  spectral  ratio  were  determined  by  assuming  the 
ship  to  be  proceeding  at  8.7  knots  into  head , 1 ong  crested  waves.  Super- 
imposed on  the  plot  are  model  test  results  from  Reference  11  for  simulated 
ship  speeds  of  25  and  30  knots  in  head  regular  waves.  Considering  the 
assumptions,  the  correspondence  is  good  between  the  model  results  and  the 
results  derived  in  part  from  the  radar  wave  spectrum.  The  magnitude  of 
differences  between  the  model  test  data  at  high  speed  and  the  ship  data 
at  relatively  low  speed  are  in  line  with  similar  results  from  model  tests 
of  a high  speed  ship  shown  in  Reference  12. 

On  the  basis  of  these  results  it  appeared  that  the  radar  data 
reduction  system  was  as  checked  out  as  it  was  ever  going  to  be  in  advance 
of  application  to  a large  number  of  intervals.  A description  of  the  final 
presentation  is  deferred  until  a later  section  of  the  present  report. 

Retrospect i ve 

A final  summary  is  perhaps  in  order  of  a few  fundamental  problems 
known  to  exist  in  the  data  and  the  procedure  and  which  have  a potential 
for  introducing  errors  of  problematical  magnitude. 
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1.  Loss  of  absolute  radar  range  --  effects  considered  of  probable 
relatively  low  magnitude  owing  to  actual  radar  aim. 

2.  Angle  measurements  influenced  by  surge  and  sway  accelera- 
tions --  effects  upon  estimates  of  true  vertical  acceleration 
considered  relatively  small  owing  to  the  nature  of  accel- 
erometers --  effects  upon  vertical  component  of  radar  range 
considered  potentially  much  more  serious  for  the  data  involv- 
ing an  accelerometer  package  not  at  the  position  of  the  radar. 

3.  Double  integration  of  very  low  frequency  components  of  accel- 
eration — the  lower  the  frequency  the  less  well  resolved  the 
the  acceleration  data  becomes,  and  the  less  reliable  the  present 
procedure. 

4.  The  assumption  that  waves  illuminated  by  the  radar  are  not  dis- 
torted by  the  ship  — there  was  no  possible  alternative  to 
this  assumption  in  the  present  work,  and  no  impression  of 
error  from  this  source  was  attempted. 

5.  The  fact  that  the  spot  on  the  free  surface  which  is  illumi- 
nated by  the  radar  moves  in  a more  or  less  sinuous  path  in  the 
horizontal  plane  due  to  the  lateral  motions  of  the  ship  — 
owing  to  item  number  2 and  the  lack  of  yaw  measurements  the 
magnitude  of  this  motion  cannot  be  fully  derived.  However, 
given  the  actual  geometry  of  the  installation,  the  amplitude 
of  such  a sinuous  path  may  be  estimated  as  about  half  a ship 
beam  at  worst,  an  order  of  magnitude  less  than  the  wave  lengths 
of  most  structural  concern.  Thus  for  long  wave  length  compo- 
nents these  effects  may  be  considered  as  of  less  concern  than 
others  just  mentioned.  In  the  range  of  wave  lengths  less  than 
ship  beam  the  distortion  introduced  into  the  encountered  short- 
wave elevations  could  be  very  large  --  to  the  point  that  the 
physical  significance  of  both  the  short  period  components  of 
the  derived  wave  elevation  and  the  corresponding  high  frequency 
spectral  densities  may  be  questioned. 

ANALYSIS  OF  THE  TUCKER  METER 

Backg round 

Since  its  introduction  twenty  years  ago,  Reference  13,  there  has 
been  considerable  effort  (Ref. 14)  expended  in  the  calibration  of  the 
"Tucker  Meter"  installations  in  the  relatively  small  (225  foot)  British 
weather  ships.  Unfortunately,  full  scale  calibration  has  apparently  not 
been  carried  out  for  larger  ships  (Ref. 15).  In  the  present  application 
to  a high  speed  ship  nearly  1000  feet  in  length  the  analyst  is  essentially 
faced  with  a background  comprising  only  the  just  mentioned  references, 
the  equipment  manuals  (Refs.  16,17)  Pierson's  analysis  (Ref. 7),  and  the 
computed  theoretical  correction  curves  attributed  to  DTNSRDC  in  Refer- 
ence 2. 
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Quoting  the  designer  of  the  Tucker  Meter  (Ref. 13): 

"In  principle,  the  instrument  measures  the  mean  height  of 
the  water  surface  relative  to  a point  on  the  ship's  hull,  and 
adds  to  this  the  height  of  the  point  relative  to  an  imaginary 
reference  surface,  thus  giving  the  height  of  the  water  surface 
independent  of  vertical  position  of  the  ship." 

The  first  height  cited  is  measured  by  pressure  gages  mounted  below  the 
waterline,  and  the  second  is  the  result  of  doubly  integrating  acceler- 
ometers mounted  adjacent  to  the  pressure  gages.  Thus  in  order  to  produce 
a practical  and  reliable  instrument  within  the  confines  of  pre-1956  tech- 
nology, some  compromises  relative  to  ideal  design  had  to  be  made.  The 
basic  description  of  the  instrument  (Ref. 13)  and  the  rationalization  of 
any  necessary  design  compromise  has  been  made  almost  entirely  under  the 
implied  condition  of  zero  ship  speed.  Indeed,  the  first  page  of  Refer- 
ence 13  includes  the  statement: 

"The  ship  must  generally  be  hove-to  when  the  recorder  is 
in  use  since  false  wave  periods  will  be  recorded  if  she  is  moving 
towards  or  away  from  the  waves;  it  has  also  been  found  by  experi- 
ment that  the  apparent  wave  heights  may  be  altered." 

Effectively,  what  was  done  in  calibrating  the  meter  for  the 
weatherships  was  to  drop  a buoy  to  obtain  a standard  wave  measurement, 
retrieve  the  buoy,  then  record  Tucker  Meter  output  while  steaming  around 
in  the  area,  and  then  finally  drop  the  buoy  again  for  a closing  wave 
standard.  Analysis  (in  R ef.lk)  involved  deriving  a directional  spectrum 
from  the  buoy  measurements,  integrating  this  along  various  ship  courses 
and  speeds  to  obtain  a predicted  encountered  wave  spectrum,  and  then  com- 
paring the  results  of  the  derivation  with  the  spectrum  of  the  Tucker  Meter 
signal.  The  Tucker  Meter  signal  was  assumed  to  be  the  result  of  cascading 
the  true  surface  elevation  spectrum  through  two  filters.  That  is,  the 
true  spectrum  Sr-(uu),  was  assumed  to  be  of  the  form: 

Scc  (uo)  = STT(u>e)/[BS(u:e)  As  ( k)  ] (50) 

where 

S—tiu  ) is  the  spectrum  of  the  Tucker  Meter  Signal 

A(k)  is  a wave  attenuation  factor  which  compensates 

for  the  fact  that  the  meter  measures  mean 
dynamic  head  at  some  depth,  whereas  the  spectrum 
of  the  surface  elevations  is  desired. 

B(uj  ) is  the  known  amplitude  response  of  the  double 
integrators  built  into  the  meter. 

Effectively,  if  the  true  spectrum  is  known  from  the  buoy  records 
the  factor  A(k)  (the  quantity  "k"  is  wave  number)  may  be  estimated  by: 
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(51) 


A2  (k) 


STT(u,e) 


SCC(u>e)B  (a>e} 

From  simple  theoretical  considerations  the  form  of  A(k)  is  expected  to  be 

-kD 


A ( k)  = e 


(52) 


where  0 is  the  effective  mean  depth  of  the  pressure  recorders  and  k 
corresponds  to  a mean  frequency.  A best  fit  of  rather  badly  scattered 
data  for  I /A  (k)  implied  D = 17.5  feet  for  the  weathership,  a value  some 
2-1/2  times  that  of  the  actual  value  of  submergence  of  the  pressure  heads. 
This  discrepancy  is  attributed  to  interference  with  waves  by  the  ship. 

The  scatter  in  derived  values  of  I/A'  (k)  ranged  up  to  a factor  of  5, 
factors  of  2 being  typical.  Use  of  the  average  correction  factor  on  large 
numbers  of  Tucker  Meter  Spectra  from  the  weatherships  is  regarded  as  being 
statistically  correct.  The  results  of  the  individual  calibrations  in 
Reference  14  imply  that  with  the  best  available  information  the  correction 
of  any  individual  Tucker  meter  spectrum  may  easily  result  in  errors  in 
spectral  area  of  a factor  of  2.  The  scatter  in  the  correction  factor 
increased  substantially  with  forward  speed,  which  no-doubt  led  to  the 
practice  of  discarding  or  viewing  with  suspicion  weathership  wave  data 
obtained  with  the  meter  if  ship  speed  was  in  excess  of  about  5 knots. 

The  top  speed  in  the  trials  analyzed  in  Reference  14  was  14  knots. 


The  computation  method  given  in  Reference  16  for  the  Tucker  Meter 
response  (or  correction  factor)  incorporates  in  an  addendum  the  results 
of  the  analysis  of  Reference  14.  in  effect,  the  wave  pressure  attenua- 
tion factor,  Eq.  (52)  is  computed  for  a depth,  D,  equal  to  2-1/2  times 
the  actual  depth  of  the  installed  pressure  sensors,  and  the  result  is 
multiplied  by  the  double  integrator  amplitude  response.  The  correction 
curves  given  in  Reference  2 appear  to  have  been  derived  in  this  manner. 


The  SL-7  Tucker  Meter  Installation 

In  the  present  case  the  Tucker  Meter  was  installed  at  Frame  119, 
the  pressure  taps  being  located  l6  feet  below  the  34  foot  waterline  port 
and  starboard.  The  beam  of  the  ship  in  way  of  the  meters  is  105.5  feet, 
and  it  may  be  noted  that  the  Tucker  Meter  installation  is  470  feet  aft 
of  the  radar  installation,  and  approximately  165  feet  aft  of  midship. 

The  prescribed  calibration  procedure  (Ref. 17)  was  carried  out 
upon  installation.  Though  this  procedure  is  highly  detailed,  one  fairly 
important  detail  is  omitted.  This  is  the  sense  of  the  signal.  As  noted 
in  Reference  5 no  sense  information  was  available  in  the  calibration  phase 
of  the  present  project.  Accordingly  it  was  necessary  to  infer  the  sense 
from  other  channels.  Of  all  the  other  channels,  that  with  the  best  reso- 
lution and  sensitivity  to  wave  position  is  the  longitudinal  vertical 
bending  stress.  A positive  stress  is  a tension  in  the  deck,  a hogging 
moment.  From  previous  data  as  well  as  the  SL-7  model  test  data  (Ref. 11) 
it  is  expected  that  a hogging  moment  would  be  in-phase  with  wave  crest 
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near  or  aft  of  midship.  Accordingly  in  the  present  case,  positive  longi- 
tudinal bending  stress  would  be  expected  to  be  roughly  in-phase  with  the 
Tucker  Meter  crest.  In  the  sample  time  histories  in  References  2 and  3 
and  those  produced  for  the  sample  run  401  previously  described,  positive 
stress  and  positive  Tucker  Meter  signals  were  roughly  out-of-phase. 

Thus  it  appeared  that  in  the  data  as  originally  calibrated,  a positive 
Tucker  signal  corresponded  to  a wave  trough,  and  that  the  sense  would 
have  to  be  reversed  in  order  to  be  consistent  with  the  crest  positive 
convention  adopted  for  the  radar.  This  is  a simple  sign  reversal  in  the 
data  processing. 

Ana!  of  Possible  Corrections 

In  the  present  data  set,  ship  speeds  are  usually  30  knots.  Wave 
directions  and  the  complexity  of  the  wave  pattern  are  known  only  to  the 
extent  of  the  few  words  in  the  log  books.  While  the  speed-length  ratio 
range  for  the  SL-7  is  about  the  same  as  the  range  tested  in  the  cali- 
bration of  the  weathersh i ps , most  of  the  present  data  involves  speeds 
corresponding  to  the  uppermost  speed-length  ratio  of  the  weathership 
trials  --  the  situation  where  the  scatter  was  worst.  Though  the  SL-7 
is  four  times  the  length  of  the  weatherships  the  pressure  gages  are 
located  only  at  twice  the  depth.  As  ships  go  the  SL-7  is  quite  a dif- 
ferent animal  than  the  weathership. 

It  may  be  noted  that  in  the  case  of  the  weatherships  the  important 

range  of  storm  wavelengths  in  the  North  Atlantic  (300  to  600+  feet)  is 
entirely  above  the  length  of  the  ship  while  in  the  case  of  the  SL-7  it 
is  entirely  below.  Thus,  for  this  important  wavelength  range  the  ship- 
wave  interaction  for  the  SL-7  would  be  expected  to  be  entirely  different 
than  that  for  the  weatherships.  Many  model  tests  in  short  oblique  waves 
have  shown  very  large  attenuation  of  incident  wave  height  along  the 
ship  --  relatively  a much  larger  attenuation  than  is  experienced  in  waves 
longer  than  the  ship. 

There  seems  no  justification  for  the  assumption  that  quantitative 
correction  factors  suitable  on  average  for  the  weatherships  would  be 
applicable  to  the  SL-7  for  the  prevalent  wave  lengths  in  the  North  Atlantic. 

For  initial  data  reduction  purposes  it  was  clear  that  at  least  the 
scalar  spectra  of  the  Tucker  Meter  signal  would  have  to  be  produced.  It 
was  evident  that  some  inspired  guessing  would  have  to  be  done  if  wave 
attenuation  corrections  of  some  sort  were  to  be  applied  to  the  data,  and 
it  was  decided  not  to  attempt  this  in  the  basic  data  reduction  procedure. 

A Study  of  Errors  Induced  by  Double  Integration 

There  are  two  aspects  to  the  correction,  Eq.  (50).  With  the  type 
of  data  in  hand  there  is  always  the  possibility  of  making  a linear  cor- 
rection dependent  only  on  encounter  frequency;  that  is,  the  correction 
for  the  double  integrator  frequency  response  was  theoretically  feasible. 
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With  the  philosophy  that  half  a correction  might  be  better  than  none, 
the  details  of  the  Tucker  Meter  were  re-examined. 


After  re-examining  the  descriptions  of  the  system  (Ref. 13,16, 

17)  it  appeared  that  the  representation  for  the  response  of  the  instru- 
ment, as  advanced  by  the  designer  and  reflected  in  Eq.  (50),  is  incon- 
sistent with  the  operations  which  are  supposed  to  be  performed  on  the 
observations.  It  is  claimed  that  the  response  is  the  product  of  a wave 
frequency  response  and  a double  integrator  frequency  response.  When 
frequency  response  is  stated  in  this  way  it  means  that  the  operations 
are  cascaded:  in  the  present  case  that  the  wave  attenuation  applies 
also  to  the  displacement;  and  that  the  double  integrator  response  applies 
also  to  the  pressure  head.  Neither  seemed  true  according  to  the  stated 
principle  of  operation  and  the  description  of  the  system. 

The  above  is  the  same  point  of  view  which  was  taken  by  Pierson 
(Ref. 7)  (though  it  was  not  so  baldly  stated  by  him),  and  in  following  up 
this  alternate  viewpoint  Pierson's  development  in  Reference  7 was  followed. 


Figure  15  indicates  the  geometry  at  the  section  where  the  meter  is 
installed.  A positive  roll  (<p)  and  a positive  vertical  displacement  Z(t) 
are  indicated.  The  two  pressure  gages  are  located  a distance  0 below  the 
nominal  load  waterline.  Associated  with  each  gage  is  an  accelerometer 
which  is  gymbal  mounted.  It  appears  from  the  descriptions,  Reference  13, 
that  the  accelerations  observed  will  be  quite  good  approximations  to  true 
vertical  accelerat ions . Accordingly,  the  quantities  actually  measured 
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Hs)  and  dynamic  (Fp,  F^)  heads  at  the  location  of  the  pressure  gages 
are  noted  in  the  figure. 


By  working  through  hand  books.  References  16,17,  it  appeared  that 
what  is  done  to  the  measured  signals  by  the  computing  circuits  in  the 
instrument  may  be  written  as  follows: 

T (t)  = 1.2[Pp-D]/2  - 1 .2A(T)*Zp(t-x)/2 

+ 1.2[Ps-D]/2  - 1 .2£(t)*Z*s  (t-x)/2  (53) 

T ( t ) is  the  indicated  wave  elevation 

■£(t)  is  the  impulse  response  of  the  double 
integration  circuits. 

The  asterisks  in  Eq . (53)  denote  the  linear  convolution  operation, 
and  the  dimensions  of  the  terms  correspond  to  pressure  head.  In  Eq.  (53) 
the  negative  sign  on  the  doubly  integrated  accelerations  comes  about 
because  the  displacement  Z(t)  is  here  defined  as  positive  downward,  while 
wave  (crest)  is  positive  upward.  The  factor  of  1.2  is  a factor  introduced 
in  the  calibration  procedure  to  compensate  in  part  for  the  apparent  atten- 
uation of  output  which  is  implied  by  the  nominal  computed  response  of  the 
i nstrument . 


where: 
and : 
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FIGURE  15  DIAGRAM  OF  THE  TUCKER  METER  INSTALLATION 
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As  may  be  noted,  the  value  of  Eq . (53)  is  zero  for  the  case  of  no 
motions  and  no  waves  since  each  pressure  is  then  just  equal  to  static 
head,  D.  In  the  same  way  if  the  integrators  were  perfect  and  good  to 
O.C.,  the  value  of  Eq.  (53)  would  also  be  zero  for  no  waves  and  the  case 
that  the  ship  was  given  a static  displacement. 


Because  of  the  accelerometer  mountings  mentioned  and  the  linearity 
of  the  double  integration,  the  two  convolution  terms  in  Eq . (53)  may  be 
replaced  by: 

-1.2*(t)*Z(0 


where: 


2 (t)  = [Z  (t)  + Zs(t)]/2 


= the  mean  vertical  acceleration  at 
the  ship  section  in  question. 


Similarly  let 

P = [Pp+Ps]/2 

= the  mean  observed  pressure 


Then 

T(t)  = 1.2P  - 1.20  -1.2Z(T)**Z*(t-T)  (54) 

Equation  (54)  represents  what  the  instrument  apparently  does.  To 
proceed  further  it  was  assumed  that  the  pressures  could  be  represented 
as  follows: 


P = F + H 
P P P 

ps  ■ Fs  * Hs  (55) 

In  Eqs.  (55)  F and  F^  represent  the  dynamic  contributions  due  to  waves 

and  the  interference  of  the  ship,  and  H and  H are  the  static  heads  from 

P s 

nominal  mean  water.  All  these  quantities  are  functions  of  time.  It  should 
be  remarked  that  the  arbitrary  division  between  dynamic  and  static  in 
Eq.  (55)  is  consistent  with  the  usual  small  amplitude  wave  theory  where 
the  total  pressure  at  depth  is  the  sum  of  a dynamic  component  due  to  waves 
and  the  hydrostatic  pressure  at  depth.  Now  computing  the  hydrostatic  heads 
from  Figure  15  and  substituting  for  the  port  and  starboard  pressures: 


P = [ P + P 1/2 
P s 

- TFp  + Fs]/2  + [Hp  + Hs]/2 

= fF  + Fs]/2  + Z(t)  +Dcoscc  (56) 


Substituting  Eq.  (56)  into  Eq.  (54),  the  expression  for  the  Tucker 
Meter  signal  becomes: 
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e(t)  = £(t)*2( t-x)  - Z(t)  (59) 

The  approximation  in  the  term  involving  D is  within  1%  for  the  range  of 
rolling  angle  involved  in  the  present  data.  In  any  event,  the  term  0^/2 
would  seldom  exceed  one  foot  in  magnitude  in  the  present  data  set. 

Because  this  term  is  quadratic,  its  contribution  to  f 1 uctuat ions  would  be 
half  a foot  at  most,  a magnitude  which  must  probably  be  considered  less 
than  the  probable  best  accuracy  of  the  system. 

If  the  ship  is  stationary,  Eqs.  (58)  and  (59)  say  that  the  "mean 
dynamic  head",  F(t),  is  conceptually  the  same  as  the  Tucker  Meter  signal, 
T(t),  since  it  should  be  remembered  that  the  factor  of  1.2  arose  in  the 
calibration  as  a gain  bias  meant  as  an  average  compensation  for  the 
response  of  the  system.  If  the  double  integration  is  perfect,  s(t)  is 
zero,  and  the  "mean  dynamic  head"  and  the  Tucker  Meter  signal  differ  con- 
ceptually only  by  a quadratic  term  in  roll  already  considered  largely 
neg 1 i g i b 1 e. 

Equation  (58)  thus  expresses  the  outcome  of  following  the  previ- 
ously mentioned  alternate  point  of  view  regarding  the  response  of  the 
instrument.  The  correction  for  double  integration  response  appears  in 
a quite  different  way.  In  the  frequency  domain,  the  scalar  spectrum  of 
mean  dynamic  head  corresponding  to  Eq.  (58)  would  very  nearly  be  the  sum 
of  the  spectrum  of  the  Tucker  Meter  output,  the  spectrum  of  e(t),  and 
the  co-spectrum  of  the  two  terms.  Presumably,  the  spectrum  of  the  mean 
dynamic  head  would  be  corrected  by  a frequency  dependent  factor  similar 
to  A (k)  of  Eq.  (50)  which  involves  only  the  attenuation  of  the  waves  and 
ship-wave  interference. 

It  appeared  worthwhile  to  investigate  at  least  approximately  what 
the  integration  error,  e(t),  Eq.  (59)  amounts  to  for  present  data.  Taking 
the  Fourier  Transform  of  Eq.  (59): 

e(u>)  = Z(u))[-I/id?]L(id)  - Z(uj) 

= Z (oj)  [ L ( uu)  -1]  (60) 

where  the  bars  denote  the  Fourier  Transform.  The  [-1/uT]  factor  is  the 
response  of  the  perfect  integrator.  The  product  of  this  factor  and  the 
acceleration  transform  is  the  displacement  transform.  The  (complex) 
function  L ( uu)  is  the  response  of  the  double  integration  as  given  in  Ref- 
erence 13.  Relative  to  true  displacement  L(uj)  is  a peculiar  high-pass 
filter  in  that  its  phases  are  leads  at  low  frequency  and  zero  at  very 
high  frequencies. 
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Equation  (60)  suggested  an  approach  to  the  evaluation  of  e(t)  using 
the  present  data.  In  the  data  reduction  for  the  radar  an  estimate  was 
produced  of  the  time  history  of  the  true  vertical  displacement  at  the  for- 
ward accelerometers.  This  is  the  only  vertical  displacement  estimate 
possible  with  the  present  data  set.  The  true  vertical  displacement  at 
the  position  of  the  Tucker  Meter  may  be  estimated  by: 

z(t)  = zradar(t)  + 470  e(t)  (6i) 

where  fl(t)  is  the  pitch  angle,  the  dimensions  are  feet,  and  the  plus  sign 
results  from  the  convention  that  displacements  are  positive  down  and  pitch 
is  positive  bow  up.  Numerically,  Eq.  (6l)  is  not  such  a good  idea  because 
pitch  is  small,  relatively  poorly  resolved,  and  the  distance  between  wave 
meters  is  enormous.  The  alternative  was  to  do  nothing. 

Now  with  an  estimate  of  true  vertical  displacement,  Z(t),  the  fast 
convolution  method  described  in  conjunction  with  the  radar  data  reduction 
may  be  used  to  estimate  e(t),  and  this  estimate  used  with  the  Tucker 
Meter  data  to  estimate  the  mean  dynamic  head,  Eq.(58). 

The  item  not  defined  in  the  data  set  is  the  function  [l((a))-1'1  in 
Eq.  (60) . All  that  is  available  as  a quantitative  description  of  L(ui)  is 
contained  in  Figures  5 and  6 of  Reference  13.  To  try  out  the  computation 
of  e(t)  some  reasonable  interpolator  function  was  required.  Arbitrarily, 
the  function  was  assumed  in  the  form  of  a second  order  high  pass  filter 
of  the  form: 


3 , 3 3 , , . 

ui  (w  -x  ) + i ax  x 
o o 

/ 3 3x3  3 3 3 

(uj^  -u)  ) + a u)q  uj 


where  u)  = a constant  cutoff  frequency 
° 

= 2tt/T 


The  two  unknown  parameters  in  Eq.  (6l)  were  selected  so  as  to  pro- 
duce a fit  to  the  amplitude  response  given  in  Figure  5 of  Reference  14 
within  about  1%,  and  a fit  to  the  given  phase  leads  within  2 . The  final 
values  used  were: 


T = 31  sec 
o 

a = 1.605 


Only  a very  approximate  check  could  be  made  of  the  fit  for  encounter 
periods  in  excess  of  50  seconds  because  of  the  nature  of  the  graphical 
presentation  of  data  in  Figure  6 of  Reference  13.  The  very  low  frequency 
behavior  of  the  interpolator  was  quite  reasonable  relative  to  this  latter 
data.  Before  applying  the  approximation  to  actual  data,  sinusoidal  dis- 
placements were  transformed  and  the  results  analyzed  to  insure  that  the 
frequency  domain  operator  [L(uj)-1],  Eq.  (60)  had  been  programmed  correctly. 
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The  procedure  was  applied  to  the  interval  previously  selected  for 
trial  purposes.  Run  1+01.  Figures  l6a  and  1 6b  are  samples  of  the  result- 
ing time  histories.  The  portion  of  data  involved  is  the  same  as  that 
shown  in  Figures  5 and  12.  At  the  top  is  the  vertical  displacement  at 
Frame  119,  Eq.  (6l).  It  may  be  noted  that  the  trace  is  noisy  due  to  the 
noise  in  pitch,  and  that  the  first  120  seconds  show  the  effect  of  tapering 
which  is  in  the  computed  displacement  at  the  radar.  The  second  frame 
shows  the  hydrostatic  head  correction,  € ( t ) , Eq.  (59).  This  trace  is  also 
tapered  as  required  by  the  fast  convolution  method.  The  third  trace  down 
is  the  Tucker  Meter  signal,  calibrated  and  with  sense  reversed  as  previ- 
ously noted.  Finally,  the  trace  at  the  bottom  is  the  "mean  dynamic  head", 
Eq.  (58).  The  term  involving  roll  was  included  in  the  computation,  but, 
as  expected,  contributed  little. 

The  surprise  in  these  results  was  how  large  the  correction  turned 
out  to  be.  The  results  for  mean  dynamic  head  are  of  the  same  magnitude 
as  those  for  the  radar  signal.  Figure  12  (the  scales  are  the  same),  whereas 
the  Tucker  signal  is  much  smaller.  What  appears  to  produce  the  large  cor- 
rection is  the  combination  of  phases.  Analysis  of  the  study  involving 
sinusoids  indicated  that  the  value  of  e(t)  for  periods  of  20  seconds  should 
be  almost  as  large  as  Z(t).  The  magnitude  of  the  correction  shrinks  to 
20%  of  displacement  for  4 sec  encounter  periods,  and  very  rapidly  to  zero 
for  shorter  periods. 

It  may  be  noted  that  the  hydrostatic  head  correction  and  the  Tucker 
meter  signal  resemble  each  other  and  are  roughly  in  phase.  Except  that 
the  various  signals  were  recorded  at  the  same  time  on  the  same  tape,  the 
two  traces  are  derived  from  independent  measurements.  For  the  long 
encounter  periods  of  the  example  the  phase  shifts  in  the  analog  double 
integrators  apparently  results  in  significant  error  in  correcting  for  the 
stat i c head . 

Under  the  stated  conditions  for  Run  401  (head  seas,  8.7  kt  ship 
speed)  the  wave  lengths  corresponding  to  the  roughly  15  to  20  second 
encounter  periods  shown  range  between  1-1/2  and  2 ship  lengths.  If  this 
was  what  really  existed  it  would  be  expected  that  the  vertical  displace- 
ment (heave)  near  midship  would  be  roughly  equal  to  the  wave  fluctuations. 
Figure  16  indicates  that  the  mean  dynamic  head  and  the  vertical  displace- 
ment magnitudes  are  roughly  equal,  and  roughly  the  same  as  the  wave  trace 
from  the  radar.  In  addition  (noting  that  displacement  in  Figure  16  is 
positive  down)  the  indicated  wave  aft  of  midship  would  be  expected  to 
slightly  lag  negative  displacement.  It  is  evident  from  the  figures  that 
this  is  true. 

Spectra  were  computed  for  the  Tucker  Meter  and  mean  dynamic  head 
data  developed  for  Run  401.  These  are  shown  in  Figure  17  in  comparison 
with  the  spectrum  from  the  radar  measurement.  As  noted,  the  significant 
mean  dynamic  head  (4  rms)  compares  closely  with  that  computed  for  the 
radar  data.  The  spectral  shape  is  different,  the  radar  spectrum  being 
broader.  Again  assuming  long  crested  head  seas,  a wave  attenuation  cor- 
rection factor  was  derived  from  Eq.  (52)  with  D = twice  the  actual  depth 
of  the  pressure  gages,  and  this  was  applied  to  the  mean  dynamic  head 
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FIGURE  17  - SPECTRA  DERIVED  FROM  "JCKER  DATA  (RUN  401) 


spectrum.  In  this  case  ths  wave  attenuation  correction  is  not  suffi- 
cient to  align  the  mean  dynamic  head  and  radar  spectra.  The  computed 
Tucker  Meter  correction  curves  from  Reference  2 are  roughly  unity  in  the 
encounter  range  of  significance  in  Figure  17  so  that  the  large  disparity 
between  the  basic  Tucker  output  and  the  other  spectra  cannot  be  resolved 
in  this  way. 

It  was  felt  that  further  results  for  "mean  dynamic  head"  might  be 
of  interest,  and  it  was  determined  to  include  this  in  the  Tucker  data 
reduction  procedure.  This  is  of  course  only  half  a correction  from 
another  point  of  view.  It  is  not  known  how  to  make  the  wave  attenua- 
tion/interference correction  to  either  the  Tucker  Meter  or  mean  dynamic 
head  data. 


INDIRECT  METHODS  OF  WAVE  MEASUREMENT 

Because  linear  systems  are  •'eversible  in  the  data  reduction  sense, 
it  is  tempting  to  attempt  an  inferential  procedure  whereby  the  wave  spec- 
trum is  estimated  from  observed  response  spectra  and  the  (assumed  known) 
linear  properties  of  the  response.  To  illustrate  the  general  idea  and 
some  of  the  more  serious  complications,  several  hypothetical  cases  may  be 
formul ated : 

CASE  1 : Ship  Speed  Zero,  Waves  Long-Crested 

In  this  case  the  response  and  wave  spectrum  are  simply 
rel ated : 

Y (to)  = 1 H (uu)  S (to) 

in  which  H(uj)  is  the  response  transfer  function 
Y(u>)  is  the  observed  response  spectrum 

and 

S(tu)  is  the  scalar)  spectrum  of  waves. 

Clearly  if  the  response  spectrum  is  observed  and  the 
response  transfer  function  is  known,  S((d)  may  be  inferred. 
The  main  complication  in  this  simple  problem  is  that  most 
ship  response  functions  look  like  low-pass  filters  with  a 
more  or  less  strong  high  frequency  attenuation.  This  means 
that  the  high  frequency  response  spectral  density  tends  to 
be  driven  down  into  the  noise  level  of  the  spectral  analysis 
method  and  the  resulting  estimate  can  be  statistically  bad. 


(64) 
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CASE  2:  Ship  Speed  Non-Zero,  Waves  Long-Crested 
In  this  case,  i n genera  1 1 

YU  ) = ' ((d)  |2S,  (ui  ) 

+ I H (u>  ) I*  S («)  ) 
s e s e 

+ 1 H,  («)  |SS3  (ui  ) 


(63) 


in  which 


U)g  is  encounter  frequency 


S.(tue)  is  the  encountered  spectrum  in  region  i 


The  encountered  wave  spectrum  is  composed  of  three  parts, 
each  of  which  has  a one-to-one  inverse  to  the  wave  fre- 
quency plane.  Thus  in  the  above  case  there  are  three 
unknown  rather  than  1.  The  solution  can  be  approached  by 
considering  3.  responses  simultaneously,  thus  giving 
3 equations  in  3 unknowns,  which  may  be  solved  j_f  the 
matrix  of  responses  is  not  ill-conditioned. 


CASE  3:  Ship  Speed  Non-Zero,  Waves  Short -Crested 


Here  it  is  instructive  to  consider  multiple  responses  and 
write  down  the  general  form  of  the  cross  and  scalar  spectra 
for  responses  j and  k: 


Yjk(u)e}  = IHjV“e)Hkl  (u,e)Sl(u>e  Xe)dXe 
+ IH]2(U)e)Hk2(u)e)S2(VXe)dXe 
+ |Hj'3(u)e)Hk3(»e)S3(u,e,Xe)dXe 


in  the  above 


V»e> 


cross  spectrum,  responses  j and  k 


Hj  i (»e> 


= response  j transfer  function  for  wave 
mapping  region  i (asterisk  denotes 
complex  conjugate) 


S.(x  ,X  ) = encountered  directional  wave  spectrum, 
1 6 6 region  i 


Xg  ~ heading  angle 


Just  as  with  the  wave  spectrum  itself,  all  the  possible 
response  spectra  and  cross  spectra  involve  the  direc- 
tional wave  spectrum  imbedded  in  i nteg ra 1 equations. 

It  is  not  immediately  apparent  if  this  representation 
can  be  inverted  without  some  gross  approximations  about 
the  relative  wave  headings.  The  St.  Den i s-P i erson  model 
is  not  a conventional  linear  system  in  the  encounter 
frequency  plane. 


Given  that  much  of  the  present  data  set  involves  high  ship  speed, 
and  that  there  are  relatively  few  channels  of  information  available,  it 
was  felt  advisable  not  to  attempt  to  plan  on  a serious  attempt  along  the 
lines  just  discussed.  For  a few  hove-to  cases  in  the  data  set  the  nominal 
ship  heading  was  head  seas,  *f>  that  there  is  the  possibility  of  trying 
Case  1 on  a few  pieces  of  data.  For  this  purpose  the  midship  stress 
spectrum  is  probably  the  most  promising  response  spectrum,  and  thus  compu- 
tation of  the  stress  spectrum  was  indicated  for  the  final  data  reduction 
procedure. 


I 
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THE  FINAL  DATA  REDUCTION  AND  PRESENTATION  PROCEDURE 


The  final  data  reduction  procedure  for  all  intervals  involved: 

a.  Four  main  computation  programs,  the  last  one  of  which 

produced  a complete  file  of  results  for  each  interval. 

b.  Two  lister  programs  to  supply  immediate  indications  of 

some  of  the  results. 

c.  One  file  consolidation  program  which  produced  one  file 

for  each  voyage  leg  containing  everything  but  the  time 

histories  of  radar  wave  and  mean  dynamic  head. 

d.  Two  programs  to  generate  the  final  graphical  presenta- 
tions for  each  interval. 

Items  b through  d amount  to  bookkeeping  operations.  The  work  is  done  in 
the  four  main  computation  programs. 

The  first  computation  program  carried  out  the  procedure  previously 
described  for  the  radar.  At  its  conclusion  the  radar  wave  spectrum  and 
the  computed  time  history  were  written  in  temporary  files  as  was  the  time 
history  of  vertical  displacement  at  the  radar. 

The  second  program  involved  reduction  of  the  Tucker  data.  Both 
the  original  data  and  the  displacement  file  produced  by  the  first  program 
were  accessed.  The  procedure  described  was  carried  out  so  that  time  his- 
tories of  mean  dynamic  head  and  the  Tucker  Meter  signal  were  available. 
These  were  spectrum  analyzed,  and  all  results  written  in  a temporary  file. 

The  third  computation  program  accessed  the  various  wave-related 
time  histories  (radar,  Tucker,  and  mean  dynamic  head)  and  performed  a 
peak-trough  analysis  on  the  middle  16-1/2  minutes  of  each.  (Because  of 
the  previously  described  tapering  both  the  radar  and  mean  dynamic  head 
data  are  not  valid  for  the  first  and  last  two  minutes  of  sample.)  The 
object  of  the  peak-trough  analysis  was  to  produce  double  amplitude 
statistics.  The  zero  crossing  convention  was  used;  that  is,  a crest  was 
defined  as  the  largest  i rstantaneous  value  in  an  excursion  above  the 
sample  mean,  a trough  was  the  smallest  instantaneous  value  in  an  excur- 
sion below  the  sample  mean.  The  double  amplitude  is  the  difference  in 
elevation  between  crest  and  succeeding  trough.  In  this  approach  small 
fluctuations  are  more  or  less  ignored  if  they  are  riding  on  top  of  large 
ones.  The  results  resemble  the  double  amplitudes  which  would  be  estimated 
by  hand  from  an  oscillograph  record  except  that  the  hand  analyst  would 
probably  visually  fair  through  superimposed  noise  whereas  the  computer 
does  not.  The  effect  is  that  while  the  computer  gets  about  the  same 
number  of  double  amplitudes  as  the  human  analyst,  the  computer's  answers 
tend  to  be  higher  if  the  records  are  noisy.  From  the  double  amplitudes 
found  the  average  of  1/3  and  1/10  highest  were  computed,  and  the  position 
in  the  sample  of  the  largest  double  amplitude  was  noted.  All  results, 
including  the  actual  double  amplitudes  were  written  in  a temporary  file. 
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The  fourth  computation  program  accessed  the  original  data  and 
performed  spectrum  analyses  upon  the  midship  vertical  bending  stress  and 
roll.  It  then  accessed  all  previously  written  temporary  files  and  pro- 
duced a new  file  containing  all  of  the  results  for  the  interval.  These 
results  included  log-book  data,  results  of  the  first  analysis  of  raw  data 
(Ref. 5),  five  spectra  along  with  all  analysis  parameters,  all  results 
from  the  peak-trough  analysis,  and  the  two  new  time  histories,  the  radar 
wave  and  the  mean  dynamic  head.  These  files  were  meant  to  be  stored  on 
magnetic  tape  for  possible  future  reference. 

The  final  presentation  of  results  for  each  interval  is  contained 
on  two  charts.  Sample  charts  for  the  example  interval  used  throughout 
this  report  are  included  as  Figures  18  and  19.  The  first  of  the  two 
charts  contains  the  scalar  spectra  and  a tabulation  of  results.  The 
second  chart  involves  sample  time  histories.  Both  are  identified  at  the 
bottom  with  the  DL  run  number,  the  voyage  number,  the  analog  tape  and 
interval  numbers,  and  the  index  number  assigned  by  Teledyne. 

Referring  to  Figure  18,  the  tabulation  at  the  left  is  intended  as 
a summary  of  the  most  significant  numbers  pertaining  to  the  interval. 

At  the  top  is  as  much  of  the  original  log-book  data  as  it  seemed  reason- 
able to  squeeze  in.  This  includes  date,  time,  position,  and  ship  speed, 
as  well  as  the  visual  estimates  of  wave  and  swell  heights  and  directions. 
Directions  are  counted  from  the  bow  to  port  or  starboard  in  degrees.  The 
"sea  state"  is  apparently  the  Beaufort  wind.  The  final  line  in  the  first 
section  of  the  tabulation  includes  comments  on  visual  weather  and,  after 
the  slash,  any  other  comment  appearing  in  the  log. 

The  second  box  in  the  tabulation  involves  midship  longitudinal 
stress  results.  Only  two  of  the  many  numbers  which  are  available  could 
be  included  as  indices.  The  first  is  the  maximum  peak  to  trough  stress 
excursion  as  obtained  in  Reference  2.  The  second  index  is  the  signifi- 
cant stress  (A  times  rms)  as  derived  from  the  area  of  the  stress  spectrum 
obtained  in  the  present  reduction. 

The  third  box  in  the  tabulation  is  a summary  of  motions.  Again 
the  "significant"  motions  (4  rms)  are  indicated.  The  value  for  roll  was 
derived  from  spectrum  area,  that  for  pitch  and  accelerations  from  the  rms 
of  the  basic  data.  (Unless  there  are  significant  linear  trends  in  the 
data  the  differences  are  slight  between  "raw"  and  "spectrum"  rms.)  The 
last  three  items  in  the  list  involve  various  stages  in  the  radar  data 
reduction.  The  first  is  the  slant  range  as  recorded.  The  "vertical 
range"  is  R (t)  of  the  radar  analysis.  This  entry  is  essentially  the 

vertical  component  of  the  range  relative  to  the  position  of  the  accel- 
erometer package.  The  number  was  derived  from  the  spectrum.  The  last 
entry  is  the  significant  displacement  at  the  radar  (significant  doubly 
integrated  acceleration).  It  too  was  derived  from  spectrum  analyses. 

In  a sense,  the  table  at  the  bottom  of  the  tabulation  contains  the 
final  numerical  answers.  Items  in  the  first  column  pertain  to  the  uncor- 
rected Tucker  Meter  signal.  The  second  column  pertains  to  the  mean  dynamic 
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head  developed  in  conjunction  with  the  analysis  of  the  Tucker  meter,  and 
the  third  column  pertains  to  wave  elevations  derived  from  the  radar  system. 
The  first  row  in  the  table  is  the  number  of  double  amplitudes  found  in 
the  middle  16-1/2  minutes  of  the  sample.  Below  this  are  noted  the  maxi- 
mum height  found  and  the  averages  of  the  1/10  and  1/3  highest  double 
amplitudes.  The  final  line  in  the  table  is  the  significant  (4  rms)  height 
derived  from  the  spectral  analyses.  Ordinarily  it  is  expected  that  the 
last  two  lines  of  the  table  will  be  about  the  same. 

At  the  right  of  Figure  18  are  plots  of  the  five  computed  spectra. 

It  was  decided  to  standardize  the  frequency  scale  from  0 to  2 rad/sec. 

In  the  great  majority  of  intervals  everything  of  interest  is  contained  in 
this  range.  In  some  intervals  one  spectrum  or  another  is  non-negligible 
beyond  2 rad/sec  but  nothing  much  has  been  seen  beyond  2.5  rad/sec  for 
any  of  the  quantities  analyzed  except  in  the  stress  spectrum  where  some- 
thing may  often  be  noticed  around  the  frequency  of  the  first  mode  of 
vertical  vibration.  The  folding  frequency  of  the  analyses  is  above 
20  rad/sec;  no  aliasing  is  expected.  Reference  5. 

The  stress  and  roll  spectra  are  plotted  together.  The  vertical 
scale  is  for  the  stress  spectrum.  The  roll  spectrum  has  been  multiplied 
by  vhe  factor  noted  in  the  legend  before  plotting.  Dimensions  of  the 
stress  spectral  density  are  (kpsi  /rad/sec)  and  those  of  the  roll  spectral 
density  are  (deg" /rad/sec)  . 

All  three  wave  related  spectra  (Tucker,  mean  dynamic  head,  and 
radar)  are  plotted  together  to  the  same  scale.  The  dimension  of  the  wave 
spectral  density  is  ( feet? /rad/sec) . In  the  wave  spectrum  plot  there  is 
a vertical  (slightly  joggled)  dashed  line.  This  line  marks  the  position 
of  the  low  frequency  cutoff,  u)^,  discussed  previously  in  conjunction  with 
double  integration  of  the  vertical  accelerations.  It  is  correct  to  inter- 
pret the  position  of  this  line  as  meaning  that  the  double  integration  has 
been  done  correctly  for  higher  frequencies,  and  incorrectly  for  lower 
f requenci es . 

There  are  several  details  about  the  spectrum  analyses  which  are 
not  documented  in  the  plots  because  they  are  constant  throughout  the  data 
reduction.  First,  the  normalization  of  the  spectra  is  such  that  the 
spectrum  area  equals  variance.  All  spectra  are  derived  from  a Fast 
Fourier  Transform  analysis  of  an  8192  point  sample.  The  fundamental 
result  is  4096  spectral  estimates  of  2 degrees  of  freedom  each.  These 
estimates  are  uniformly  spaced  in  frequency  at  a delta-frequency  of 
0.00511  rad/sec.  In  order  to  improve  statistical  reliability  the  basic 
spectral  estimates  were  averaged  in  blocks  of  20  estimates  at  intervals 
of  10  estimates.  The  resulting  averages  are  thus  equi-spaced  on  the 
frequency  scale  at  intervals  of  Au>  = 0.0511  rad/sec.  This  also  means 
that  adjacent  spectral  estimates  as  shown  in  the  plot  are  not  quite  inde- 
dendent  --  to  about  the  same  degree  as  spectral  estimates  from  the  older 
autocorrelation  methods  are  not  independent. 
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As  a result  of  the  averaging,  each  spectral  estimate  has  AO  degrees 
of  freedom  associated  with  it.  Accordingly,  the  307  confidence  bounds  on 
the  spectra  shown  in  the  charts  may  be  formed  by  multiplying  the  values 
given  by  0.72  and  1.51.  Had  the  process  sampled  continued  indefinitely 
and  a large  number  of  20.5  minute  samples  been  obtained  and  analyzed, 
nine  cut  of  ten  of  these  new  estimates  of  spectral  density  would  be 

expected  to  lie  within  the  bounds  so  constructed.  The  practical  implica- 

tion is  simply  that  the  influence  of  sampling  variability  upon  the  given 
numerical  results  is  roughly  the  same  as  that  associated  with  the  result 
of  most  other  full  scale  wave  measurement  exercises. 

The  last  detail  of  the  spectrum  analysis  is  the  "total  degrees  of 

freedom"  Reference  10.  This  number  is  included  in  parentheses  at  the  end 

of  each  line  of  legend  because  it  depends  upon  the  shape  of  each  indi- 
vidual spectrum.  It  is  an  estimate  of  the  proper  number  of  degrees  of 
freedom  to  use  in  constructing  confider.ee  bounds  on  the  sample  variance. 

If  each  of  the  numbers  in  the  present  8192  point  time  histories  had  been 
picked  randomly  the  "total  degrees  of  freedom"  would  be  8191.  This  is 
not  the  case  — adjacent  members  of  all  the  present  time  series  are  highly 
correlated  so  that  the  equivalent  "random"  sample  size  is  much  smaller. 

In  the  present  data  set  the  "total  degrees  of  freedom"  (TDF)  is  expected 
to  vary  between  60  and  600.  Approximate  90%  confidence  bounds  on  the 
variances  assuming  a Normal  zero  mean  process,  may  be  constructed  by 
multiplying  the  estimate  by  two  factors  derived  from  the  percentage  points 
of  the  Chi-square  distribution.  Examples  of  the  values  of  these  factors 
are  given  as  follows: 

TDF  Factor  for  Factor  for 

Lower  Bound High  Bound 


60 

.72 

1 .32 

120 

.80 

1.27 

200 

.84 

1.17 

400 

.89 

1.12 

600 

.91 

1.10 

These  are  factors  for  the  variances.  The  square  root  applies  to  the  rms 
values  so  that  very  roughly  the  90%  confidence  bounds  on  rms  range  from 
the  sample  rms  ± 15%  for  TDF  = 60  to  the  sample  rms  ± 5%  for  TDF  = 600. 

The  practical  implications  of  these  results  are  quite  similar  to  those 
mentioned  in  connection  with  the  confidence  bounds  on  the  spectra.  There 
is  only  so  much  "precision"  obtainable  from  one  20  minute  sample  of  wave 
elevation  --  that  which  was  attained  in  the  present  work  appears  compar- 
able to  that  achieved  in  the  past  in  similar  studies.  With  respect  to 
comparisons  between  wave  meters  or  between  data  ard  predictions  of  rms 
ship  responses  there  can  be  little  just ' f icat ion  to  a concern  about 
differences  of  5 to  15%  magnitude. 

Figure  19,  the  sample  time  histories,  needs  little  explanation, 
except  perhaps  to  say  that  the  duration  of  the  sample  shown  (8-1/2  minutes) 
was  a compromise  between  a desire  to  display  as  much  of  the  16-1/2  minutes 


66 


r 


of  derived  wave  time  histories  as  was  possible  in  one  page;  and  the  desire 
to  spread  the  time  scale  out  so  that  individual  fluctuations  were  visible 
for  intervals  involving  high  ship  speed  in  head  seas.  To  produce  the 
charts  an  8-1/2  minute  portion  of  the  available  16-1/2  minutes  of  sample 
was  chosen  such  that  the  largest  radar  wave  double  amplitude  is  shown  — 
as  well  as  (if  possible)  the  largest  mean  dynamic  head  double  amplitude. 

It  may  be  fairly  asked  why  the  effort  in  producing  plotted  time 
histories  for  each  interval  was  considered  worthwhile.  The  answer  to  the 
question  is  fairly  simple.  While  the  present  data  in  its  original  analog 
form  has  been  scanned  systematically  by  eye  (Refs.  2,3,5),  the  process 
involved  oscillograph  records  with  a time  scale  of  about  15  mi nutes  to 
the  inch.  At  this  time  compression  only  a gross  idea  of  what  was  happen- 
ing can  be  formed,  no  detailed  assessment  of  the  bel i evab i 1 i ty  of  the 
data  can  be  made,  and,  most  importantly,  the  odd  malfunction  which  is 
enough  to  upset  the  spectrum  estimates  or  the  statistics  may  often  go 
unnoticed.  This  last  is  considered  most  important  in  the  radar  data. 

It  was  pointed  out  earlier  and  in  Reference  5 that  an  attempt  was  made  to 
weed  out  intervals  where  the  radar  had  evidently  lost  signal  and  re- 
established a new  reference  range.  In  this  process  only  the  most  obvious 
instances  could  be  identifiedj.no  guarantees  could  be  made  that  all 
instances  of  moderate  or  small  magnitude  had  been  eliminated. 
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